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Theoretical Pressure Distribution at Zero Lift 
at Supersonic Speeds for Slender Delta Wings 
having Fuselages of Circular Cross Section 


R. C. LOCK, M.A., Ph.D. 


(Aerodynamics Division, National Physical Laboratory) 


SUMMARY: A method is developed for calculating the pressure distribution at 
zero lift for a class of slender configurations consisting of a wing of delta type 
plan form combined with a body of revolution of the same overall length. 
Slender-body and slender-wing theory are combined to take into account the 
interference between the wing and the fuselage. An expression is found for the 
pressure at a general point on the wing or fuselage, although some numerical 
integration is still needed. This expression is simplified to obtain algebraic 
formulae for the pressure at the wing-fuselage junction and along the top of the 
fuselage. A further simplification is possible when the fuselage radius is small 
compared with the local wing semi-span. A particular example is considered 
consisting of a “Newby” wing with a parabolic fuselage; pressure distributions 
and drag curves are calculated. 


1. Introduction 


The present programme of work in the United Kingdom on slender wings for 
supersonic aircraft began in the expectation that the final configuration would be of 
the “all wing” type, without a discrete fuselage, and accordingly most of the 
experimental and theoretical work that has been done up to the present has con- 
centrated on this type of shape. More recent design studies, however, have suggested 
that such configurations may not in fact be the most suitable, and have tended to 
return to more conventional shapes in which a thin slender wing is mounted on a 
fuselage of circular cross section. Some experimental work on these designs is 
proceeding, but there does not yet seem to be available any corresponding theoretical 
treatment. The present paper is an attempt to fill part of this gap, at least as regards 
the calculation of the pressure distribution for symmetrical uncambered wing- 
fuselage combinations at zero lift. 


There is already a considerable body of experimental evidence that, for the 
all-wing configurations mentioned previously, an adequate treatment is provided by 
the so-called “slender-wing” theory—that is, slender-body theory into which the 
further simplification has been introduced that the boundary condition may be 
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satisfied in the mean plane of the wing; so that the flow may be represented, as is 
usual in the linearised theory of thin wings, by a source distribution in this plane of 
strength proportional to the local slope. Good agreement between this theory and 
experimentally measured pressure distributions has been obtained in general, in spite 
of the fact that many of the shapes tested appear to be far from “thin” in the usual 
sense, particularly over their forward portions. But when such wings are mounted 
on a fuselage of circular cross section, the theory in its present form is clearly 
inapplicable to the prediction of detailed pressure distributions without modifica- 
tion. Fortunately, the problem can be solved by the introduction of image sources 
within the fuselage. as has been done by Lighthill®, Byrd®’ and Brebner“’. These 
authors all restricted themselves to cases where the fuselage radius R, is constant; 
when this is allowed to vary, a further complication appears because the flow field 
due to the fuselage interferes with the boundary condition on the wing, so that an 
additional source distribution in the plane of the wing must be introduced to correct 
this effect. 


There is also a possibility of obtaining improvements in performance by the use 
of wings having a blunt trailing edge. This refinement does not introduce any 
additional difficulties to the problem of calculating the pressure distribution, but 
does require a modification to the formula for the wave drag. This modification has 
been considered in Ref. 4, so that it is now possible to obtain both the detailed 
pressure distribution and overall wave drag for the most general wing-fuselage 
combination of this type. The algebraic complications tend to become excessive, 
however, except for the simplest shapes—one of which is considered in detail in 
Section 6—so that for future use a more purely numerical procedure might be 
desirable. 


NOTATION 
All co-ordinates and other lengths are measured in terms of the overall length 
of the wing-fuselage combination. 


x,y,z Cartesian co-ordinates with the origin at the apex of the wing; the 
x-axis is in the direction of the undisturbed flow and the z-axis is 
normal to the mean plane of the wing 


r= /(y? +2’) 
alternative to y 
Yio", =y/S, 1/S respectively 
C=y+iz 
Zw _ half the local wing thickness 
z, value of Zw at the centre section (y =0) 
s overall wing semi-span at the trailing edge 
local wing semi-span 
R(x) fuselage radius 
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R,=R/oe 
Ry=R (1) fuselage radius at base 


S (x) cross-sectional area of combination 
(x)=7R?’, cross-sectional area of fuselage 
Sw cross-sectional area of exposed wing (see Fig. 1) 
Swe cfoss-sectional area of gross wing 
Swe cross-sectional area of part of wing blanketed by fuselage (see Fig. 1) 
Swe* see equation (33) and Fig. 1 
A (x, y)=(0/ 0x) Zw (x, y), local wing slope 


(") (see equation (6)) 


A* (x, y)=A+ RR’ 
Oy 


(x,)dn+RR = (see equation (25)) 


(x, y)= 


0 


M Mach number 
B= /(M?-1) 
U_ velocity of free stream 
Pp) density of free stream 
q=4p,U’ 
p pressure 
p. =pressure of free stream 
C,=(p-—p.)/q. pressure coefficient 
Cyr pressure coefficient due to the fuselage alone 
Cyw pressure coefficient due to wing and interference 
C,w pressure coefficient at root section of gross wing alone 


C,; interference pressure coefficient at wing-fuselage junction (equation 
(64)) 


AC, difference in pressure coefficient between the top of the fuselage and 
the wing junction 


perturbation velocity potential 
%,.%, see equations (1) and (2) 

D_ wave drag (for unit length) 

V =(overall volume) /(length)* 


V,=(volume of fuselage)/(length)’ 
Vw =(volume of exposed wing)/(length)’ 
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K, 128) D/(qV*) (see equation (82)) 
F(x) see equation (18) 
G(x) see equation (29) 
H(x) see equation (54) 
7 root thickness/chord ratio of gross wing 
¢ fineness ratio of fuselage 
a_ length of equivalent pointed fuselage (see equation (69)) and Fig. 1) 
A=/(sa) 


2. General Theory for Pressure Distributions 


A typical wing-fuselage combination is shown in Fig. 1, which also shows the 
system of axes and notation used. The fuselage apex is at the origin, O, and the 
wing has a straight unswept trailing edge at x= 1, z=0. The wing leading edge has 
the equation y=o (x), z=0, and o(1)=s is the ratio of the wing half-span to its 
length. The fuselage has the equation r= R (x) and is assumed to be of the same 
length as the wing. The area of a cross section of the fuselage by a plane 
x=constant is denoted by Sy (x)=7R*, and the cross-sectional area of the exposed 
wing is Sw (x) so that the total area is S (x)= S»+ Sw. 


(a) PLAN FORM OF WING AND FUSELAGE b) TYPICAL CROSS SECTION 


ZZ 


(C) DEFINITION OF CROSS-SECTIONAL AREAS 


FiGure 1. 
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The perturbation velocity potential Uo can, according to slender-body theory, 
be written in the form 


(x,y, Z)+0,(X),. ; (1) 
wher (x): 5” (log dé (2) 
0 
B= /(M?-1) 


and ¢, is a harmonic function of vy and z in the cross-flow planes x=constant, which 
has to satisfy the boundary conditions 


dR 

(3) 
on the fuselage r= R (x) and 

09, OZw 

t = ty) . . (4) 


on the wing, where Zw (x,y) is the wing thickness distribution. It is normally 
sufficiently accurate to satisfy the wing boundary condition (4) in the mean plane 
z=0, but it will appear later that, in order to deal correctly with the interference 
effects, it is necessary to consider the flow induced by the fuselage at the actual 
surface of the wing (z= Zw). 


These boundary conditions are to be satisfied by a suitable distribution of 
sources in the wing plane and at the centre of the fuselage. This can best be done by 
starting with the source systems appropriate to the wing and fuselage taken 
separately and then adding additional sources to allow for their mutual interference. 
Thus the exposed wing is represented initially by a source distribution of strengtht 
2A (x, y) in the plane z=0(R < |y| <@), where 


(x, y) (5) 
Ox 


A (x, y)= 
is the local streamwise slope of the wing surface; the fuselage is similarly represented 
by a line source along the axis Ox of strength 2*RR’ (x)= Sy’ (x). 


We must now consider the additional normal component of velocity induced on 
the fuselage by the wing source distribution and on the wing by the fuselage source 
distribution. In the first place the velocity normal to the fuselage due to a wing 
source at (x, y,0) can be counteracted by an image system within the fuselage con- 
sisting of an equal source at the image point (x, R*/y, 0), together with an equal sink 
at the centre (x, 0, 0). Secondly, the source 2*RR’ representing the fuselage 


+Taking the potential of a unit source to be (1/27) log r. 
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produces a radial velocity RR’/r, and this has a component normal to the wing 
surface which is approximately equal to 


y \y oy oy \ 
To correct for this effect we must therefore introduce an equal and opposite source 
distribution in the wing plane, thereby modifying the original strength A, defined in 
equation (5), to 


A* (x, y) =A (x, y) + RR’ : (6) 

This modified source distribution A*, in the plane z=0, together with its image 

system in the fuselage, represents correctly the flow due to the exposed wing and the 

mutual interference with the fuselage. Taking into account both halves of the wing 

(y > 0 and y < 0), we obtain the following expression for the harmonic part of the 

velocity potential : — 


; dy, 


2= (7) 


et) log r+ : A* (x, n) log 


R 


where (= y+ iz and A* is defined by equation (6). 


It is interesting to note that the expression (7) for %, corresponds to a total 
source strength Q (x) in the transverse plane, where 


(since the total strength of the image 
Q (x)= Sr [ (x, 9) dn sources is zero) 
R 


Sy (x) 4 
R 


OZw 
= (x, ») dn -4R’zw (x, R) (from equation (6)) 


if 
=S, (X)+4 | zw (x, dn 
R 


= Sp’ (x) + Sw’ (x), 


since Sw (x)- 4| Zw (x, ) dy, 
R 


provided that the wing is thin relative to the fuselage. This is, of course, as it 
should be but, if the modification to A given by equation (6) had not been used, the 
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total source strength Q(x) would have been in error (by the amount 4R’zw (x, R)). 
The reason for this discrepancy can be seen more clearly by considering the harmonic 
flow in a typical transverse plane. The fuselage source (of strength 2xRR’) produces 
an outward radial velocity R’ at the circle r= R representing the fuselage, as required 
by the boundary condition there; but the total flux across that part of the circle 
external to the wing is not 2=RR’, as it would be for the fuselage alone, but 2*RR’ 
—4R’zw (x, R). This deficiency 4R’zw (x, R) can be thought of as flowing across the 
two small arcs of the circle blanketed by the wing, and thereby disturbing the wing 
boundary condition, as already described. 


It might be imagined that the neglect of the effect just described would lead only 
to a very small error in the resultant pressure calculations. It may be shown, how- 
ever, that this is not the case, and that in fact the error would be of the same order of 
magnitude as the other interference effects. 


The form taken by equation (7) when R is constant has been obtained previously 
by Lighthill®’, Byrd® and Brebner®’. 


It is convenient to rewrite equation (1) in the form 
+ Ow, 
where 4; is the potential due to the fuselage alone, 
x 


(x) log 48r— Sp” (€) log (x — dé (8) 


Pp 


0 
and $w is the additional potential due to the wing, allowing for interference effects, 


x 


ow= Sw” (€) log 2(x-8) d&é+ A* (x, log | | dyn. (9) 
0 R 
The pressure coefficient C, is given in the usual way by 
c,=-2 5° (5°) - (3 


Before proceeding further, it is desirable to enquire to what extent we should take 
into account the quadratic terms (09¢/dy)*+(éo/0z)? in equation (10); these are 
normally included in slender-body theory but are neglected in “thin-wing” theory. 
To do this it is necessary to examine the orders of magnitude of the various terms 
involved. Suppose that the wing thickness ratio and slopes are of order + and that 
the fuselage thickness ratio and slopes are of order ¢; the ratio of wing semi-span 
to length is, of course, s. Then the relevant orders of magnitude are given by Table I. 
The various pressure terms are thus of the orders of magnitude given in Table IT. 
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TABLE I 
Quantity Order of magnitude 
Sy TS 
Sy 
ox, 4 log 
F Ox 
Oop 0%, 
=—_> > ¢ (on and near the fuselage) 
or oy @ 
owe = ts logs 
cx 
Obw 
dy’ ez 
TABLE II 
Wing Linear (Ow / Ox) rs log s 
Quadratic +(0¢/0z) r2 
Fuselage Linear 0x) 2" 
Quadratic (0¢,/ dr? 
Interference rz logs 


Provided that 7s and <* are of comparable size, which will be the case when the 
wings and fuselage have comparable cross-sectional areas, then the pressure con- 
tributions due to the fuselage (both linear and quadratic terms) and the linear wing 
term will be of equal importance, at least on and near the fuselage. On the other 
hand, the restriction to “thin-wing” theory for the wing means that the quadratic 
wing terms must be neglected, so that = should be small compared with s log s and < 
if the theory is to be strictly applicable. Physically this means that the wing should 
be “thin” compared to the fuselage, even in transverse cross section; but the results 
of experiments on all-wing configurations suggest that this may not prove a serious 
restriction. 


The argument just given shows that it is legitimate to express the pressure 
coefficient as the sum of a contribution C,,» due to the fuselage alone, which should 
be obtained by full slender-body theory including the quadratic term+, together with 
a linear term C,w due to the wings and interference effects. 


2 
Dy Opp 
with Cyp=—25% ( 


‘\ 2 
(x) log (48r) — Sp” (0) log x — | Sp” (€) log (x ) 


(12) 


tIt might in fact be desirable to use more accurate methods for Cyp: particularly off the 
fuselage, since slender-body theory is only strictly valid in the immediate neighbourhood of 
the body concerned. 
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Ox 
(x) log (48) Sw” (0) log x - [ Sw” (2) log 
J 
(13) 
&) 
“2 2\(*2,2_ 
where | AX (x, 7) log 
R&® 
here (=y+iz as before. 
Thus on the wing, at the point (x, y, 0), 
ol (vy? — — R*) A® (x, n) 
(x, n) log dyn—4R*R (x) | (vy? dn+ 
R R 
272 _ 2_ 
+o” (x) A* (x, 0) log _ aR’ (x) (x, R) 5) 
yo" 


Equations (11) to (15) provide a general expression for the pressure due to a wing- 
fuselage combination of the type under consideration. All the terms except the last 
(equation (15)) are already familiar in this subject. In particular, Lord and Brebner™ 
have shown that, for area distributions which can be conveniently expressed as power 
series in x, so that 


N 

n=0 

then (£) log (x — £) dE= {S” (x) —S” (0)} logx—F (x), . (7) 
N 

where F(x)= o,7,x"*? 
n=0 
1 ) 


In fact, it may easily be seen that 
: 


this may be used to define F (x) for more general types of area distribution. 
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Equations (10) and (11) may therefore be written in the form 


l B ol 
and — [Ss (x) log Fw (x)+2 if] (20) 


The only additional complications introduced by the present problem are thus 
concentrated in the term @//0x (equation (15)); but they turn out to be very consider- 
able, except in the simplest examples. 


3. Pressure Distributions at the Wing-Fuselage Junction 


3.1. GENERAL 


In cases where the wing slope A (x, 1) can be expressed as a polynomial in 7 
(so that A* consists of a polynomial together with a term in »~?), it is possible to 
express 0//0x (equation (15)) explicitly in terms of elementary functions; but for a 
general position on the wing the results are so complicated that it would seem more 
profitable to use purely numerical methods of evaluation. It is not proposed to 
pursue the matter further in the present paper; but fortunately it is possible to 
achieve a considerable simplification as regards the pressure distribution along the 
wing-fuselage junction, which is probably the most important region from a practical 
point of view. 


If we take the limit of equation (15) as ( —> R, we find the following expression 
for 0] / 0x in the junction 


al 


| Ax (x, log dn + 20’ (x) A* (x, log 
R 


—2R’ (x) (x, n) log dn + A* (x, log +A* (x, R) log 4R | 
R 

(21) 


This expression can be simplified by separating out the terms involving 
logarithms, which can be done in the following way. Consider first the integral 


oA* n+R 


a 
Now | (1) log (1—R) dn = | (x, 1) —A* (x, R| log (1 — R) dn, 
R R 
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and, integrating by parts, this becomes 


{A*® (x, 7) (x, R)} log (o — R)- ] 


R 


A* (x, n) —A* (x, R) dn; 


(22) 


the integrand is now non-singular at »=R and, if z(x,) can be expressed as a 
polynomial in », so also can the integral. Similarly, 


OA 


~ 


* 
(x, log + R) dyn = {A* (x, 7) —A* (x, — R)} log (a + R)- 


(x, R)—A* (x, — R)} log (23) 
R 


For this equation to have a meaning it is necessary that the analytic continuation of 
A* (x, n) for negative values of » should be known; there is no difficulty here if z is 
given as a polynomial in ». Note that for wings where z is not an analytic function 
of » over the whole range (—», 7) (in particular when the cross section has a 
discontinuity in slope at 7=0), A* (x, —») is not simply the value of A* at the point 
(x, —n) on the actual wing. For example, if the wing has rhombic cross sections, so 
that z=a+ b}n|, then 


thus A* (x, — R) =a’-b’R- ‘ 


while the value of A* at the point (x, — R) of the wing is the same as at (x, + R), 
namely 


Combining these two results, (22) and (23), we find that 


| an (x, n) log dn + A* (x, log (FR +A* (x, R) log 4R 
R 


=A* (R) log (o — R)—A* (— R) log (o + R)+ A* (R) log 2+ A* (— R) log 2R + 


R 
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It is convenient to express A* as the sum of an even and an odd function of »; 


thus (x, 1) =Ag*® (x, 1) + (X, 9), 
where =F [A* (x, 1) +A* (x, —)] 


Ao® =4 [A* (x, n)—A* (x, — n)]. 


(In the example quoted, =a’ —aRR’/»;? and =b’s). 


The right-hand side of equation (24) then takes the form 


+ foe R) log 


7 
of R{A* (x, ) —Ag® (x, R)} — (x, R) dn, (24a) 
1? —R? 
R 
In order to apply the same process to the integral 
) 
Ox (x, n) log (= dn, 
R 
it is first necessary to define the function 
RR’ 
A* = ) + zw 
(x, 7) | (x, y) dy (25) 
0 
so that 
On 


where A,* and A,* are even and odd functions of respectively; clearly 


A.* 
and =A,*. 
on 
106 
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Applying the technique described, we find that 
oA* n? — R?\ o? — R? 
< / * 
| ax (x, n) log | dy + r’ (x) A* (x, log 


oA o - 
* 
(x, 7) + o’A* (x, log ( 


ax (x, R) log ( R (R) log ( 
OA* Ao* 
(x, = (x, R)! (x, R) 
| ox Ox | 0 oA yan 
(n? — R®) 
R R 
(27) 
1 
, Substituting from equations (24a) and (27) into equation (21), and using the relations 
d cA* 
dx (x, (x, (x, oc) 
(x, R)| (x, R)+ R’Ao* (x, R) 
dx Ox 
(x, R) (x, R) R’A,* (x, R), 
dx Ox 
we find that 
al; -2 * — R? 
| A® (x, log | 2k Ag® (x, R) log (7 R - 
8 
2k E (x.R)| log + 4G (x), (28) 
where 
oA 
ax (x, (x, R)} RR’ 
(x,R) —R’do* (x, R) 
x 
G (x)= 8 (n? — dn+ 
oA* dn 
+4 | (29) 
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(Note that in cases where z (x, 9) is not a purely even function of », so that Ay* + 0, 
there may be an additional term in log (c/ R), due to the last term in equation (29); 
it does not seem desirable to separate this, since no additional simplification is 
achieved in the final expression for the pressure.) 


If we define a gross wing cross-sectional area 


Swa (x)=4 [ (x, 9) dn, 30) 
0 

and write Swr — Sw. ° (31) 

so that Swr (x)=4 | Zw (x, dy @ 


0 


is the wing area blanketed by the fuselage, then it is easy to show that 


A* (x, 7) =} (x) 
and A* (x, R)=4 Swe’ (x). 


In cases where z is not a purely even function of », it is necessary to define an 
additional area 


R 
Swe® =4 | {Zw (x, —1)—Zw (x, )} dn; (33) 
then A "hel (x, R) Swr*’ (x). 


The physical definition of these areas is shown in Fig. 1. 


In terms of the areas Sw, Swr and Swp*, the final expression for the pressure at 
the wing-fuselage junction is found from equations (19) and (20) to be 


where Cyr|s= & (x) log — . (35) 
Bo R’*\|\ + R)?| 
d Ow { + 
an J E (2) log (1 =} + Sw (x) log | 4cR | 
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Here Fy (x) and Fy (x) are defined by equation (18), using the appropriate area 
distribution, and G (x) is defined by equation (29). 


When z can be expressed as a polynomial in »* only, so that A* is a purely even 
function and A* is purely odd, equation (36) simplifies further to 


Cow J | Sw” (x) log |! t Swe (x) log | 4cR + Fy (x) + G (x) 
(36a) 
oA* oA* 
(x,)-R (x, R)+ RR’ {A* (x, n) —A* (x, R)} 
with G(x)= —8 — dn+ 
R 
oA* dn 
ax (x, (29a) 
R 
N N 
If zw (x,»)= a,(x) ny", then it can be shown that G(x)= 2G, (4), 
n=0 0 
where 
n+1 a, dx 


(o"-'— R" (aRR’) (n~l), (7) 


Qn Rati l 


1 dy if n is even 


R 


dn if n is odd. 


The first three functions G,, G, and G, are given by 


” R d ata 
G, = ~—4(¢-R) E R ) (writing a@,=Z,), 
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| 


a LOCE 


d o 
—g." — R? 
G,= -a,” (¢ R*)+4 (a, RR) log 
{d 


2 
and (1+ +703) a.” 
9 o o? \dx 


(a.RR’) + 2a/RR’| 


In all the preceding work, no assumption has been made about the wing cross- 
sectional shapes except that z, should be expressible in the form of a finite poly- 
nomial = a, (x) 7". Thus equation (36), used in conjunction with standard techniques 
of polynomial fitting, could form the basis of a numerical method for calculating the 
junction pressure distribution for a wing of arbitrary shape. It is not proposed to 
proceed further in this direction in the present paper; but in the two following 
sections explicit formulae will be given for two types of wings having “similar” 
cross sections—rhombic and parabolic respectively. 


3.2. RHOMBIC CROSS SECTIONS 


The wing thickness distribution has the equation 


zw (x, (1-2), 

Sw=2z0 (1 =) (41) 

Sur (1-4 (42) 

and Swr* = = (43) 
Substituting in equations (36) and (37) (with a, =z, and ad, = —z,/o), and combining 


the term in log (/ R) in G, with the other logarithmic terms, we obtain 


| 4R, | 


(1 -R,»)! +4(z,RY’ log 


Ss ” (x) log 


Bo 
2x 


l 
+2 dx | top + Fw + 


+ (0? R’) +4R(1-R) (44) 
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in this equation, R, (x)=(R/o); Sw (x) is given by equation (41), and Fy (x) is 
derived from it by means of equation (17) or (18). 


3.3. PARABOLIC CROSS SECTIONS 


The wing thickness distribution has the equation 


so that Ning? (2:) — 
y 
, 2,RR z,\’_,. 2, RR’ 
4 
Swe = 4z.R —4R P ‘ (49) 
and Swr* = 0. 


In this case we obtain 


2 2 
: (x) log (1 - R,2)| +4 d fa R,)’| 


{z.R (1 —4R,’)} log | | 


+ Fy (x)+G (50) 


where 

4. Pressure Distributions along the Top of the Fuselage 


It is also possible to obtain relatively simple expressions for the pressure 
distribution along the top of the fuselage (v= 0, z= R). 


Writing (=iR in equation (14), and differentiating the result with respect to x, 
we obtain equation (52) which followst 


+See the footnote on p. 112. 
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LOCE 


(x, ") dy 


R R 


+ 20’A* (x. log 


—2R’A* (x, R)log2R.  (52)' 


Making use of the function A*(x, 1) defined by equation (25), such that A* =dA*/dn, 
this may be written 


tag (22) | 
A* A* 2R+ 
72 gx ( 2“ !Ax%(x, R)! log 2R 


dx | | 
¢ RR’A* (x, n)— 3 (x, 
R R 


If we compare equation (53), which gives the value of 0//0x at the top of the 
fuselage, with equation (28), which gives the value at the wing-fuselage junction, we 
notice that many of the terms are similar; and all the other terms in the equations 
for the pressure ((19) and (20)) are identical. It is therefore possible, by subtracting 
equation (28) from equation (53), to obtain the following expression for the 
difference in pressure between the top of the fuselage and the junction: 


AC, =C, (top) —C, (junction) 


A* 
| RR’X* (x. 1) ~ S— (x0) 
a Cx 
where H (x)-8 


R 
A 


) 01 
"Vax 


(x.R)) + RR’{A*(x,7) (x, R)} - (x, R)—R’A,* (x. R) 
dy. 


R? 


The first integral in H (x) involves. in general, terms in log [(o? + R*)/2R?] and in 
tan~' [(7— R)/(a + R)]; the second integral has already appeared in the function 
G (x) (equation (29)). 


+The same result may also be obtained by replacing y? by —R? in equation (15); this is per- 

missible because { occurs in equation (14) only as ‘2, which is real in both cases. Such a 
procedure is not possible for a general position on the fuselage, when the modulus in equation 
(14) must be evaluated first and then differentiated with respect to x. 
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In the case of wings with rhombic cross sections (cf. Section 3.2), it is found that 


H (x)=8 (Rz,” (7— 4.2 | og 


o? + 


2 (Rz,y" tog (7 2 (2) "| tog 


| 
+ 4(Rz,” —z,R”) tan-' 


For wings with parabolic cross sections (cf. Section 3.3), the corresponding 
expressions are 


” 


+ tan! (753) - 


5. Approximations Valid when R/c is Small 


In many cases of practical interest the fuselage radius will be fairly small 
compared with the local semi-span over much of the length; we therefore proceed in 
this section to examine what simplifications can be achieved in the theory of 
Sections 3 and 4 by assuming that R, = R/c is small. 


Considering first the pressure in the wing-fuselage junction (Section 3), we notice 
that the principal term (the first) in C,w!; is of order (z,0)” log [8o/(2x)]. If we 
neglect relative to this all terms of order R,” log R, or lower, we obtain 


ex = | Se" (x) log (57) + Swe” (x) log (45 + Fe | (57) 


since Swp” is of order R, and Swr*” of order R,’, relative to Sw”. 
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If we now examine the series (37) for G (x), we see that 


G, = d (1+0(R,?)] 


dx 
(58) 
G,=G,{1+ O(R,? log R,)} 
and (n> 2) 
A m+. 
where G,= - ip 


is the value of G, when R=0. Now, in the notation of equations (30) and (31), 


Sw = Swe Swr. 


x 


and if we define Fwr= [ 


Sw (x)- Swr” 


59 
x-€ 


dé, (60) 


then clearly Fw=Fwoe—F wy. . (61) 


If we consider the gross wing alone (continued through the body), and denote 
the value of the pressure coefficient at the centre section (v=0) of this wing by 
Cyw, then 


Cow? = (x) log (=) + Feo (| (62) 


and we can therefore write equation (57), making use of equations (58) and (61), 
in the form 


‘ 


where Cu = [Swe (x) log ( ) + R | ‘ . (64) 


To the same order of approximation Sw», which is the cross-sectional area of that 
part of the gross wing blanketed by the fuselage, may be taken to be 4z,R, so that 
the interference pressure (equation (64)) may also be written in the form 


28R 
Cyr == (x) log ( x )| (64a) 
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The final approximate expression for the pressure at the wing-fuselage junction 
is thus 


with the three terms on the right-hand side given by equations (35), (62) and (64) 
respectively. It will be noticed that the equation for the interference pressure has a 
certain similarity to the equation for the pressure due to a body of revolution whose 
area distribution is equal to Swy; this might perhaps have been expected, but the two 
expressions are not identical and the similarity can not therefore be exploited further. 


A similar simplification can be achieved by making the same approximation 
with regard to the pressure distribution along the top of the fuselage (Section 4). 
Neglecting terms of order R* log R or smaller in equations (55) and (56), we obtain 
in both cases 


AC, =C, (top) —C, (junction) ~ Sep” log 2-2 (Rz,”-z,R”). . (66) 


with Swe ~ 4Rz, as before. 
This result is also true, in fact, for wings of general cross-sectional shape. 


The approximate equations (64) and (66) are valid for a wing of quite general 
thickness distribution, and thus provide a rapid method of estimating the pressure 
at the top and junction of any wing-fuselage combination of this type once the root 
pressure distribution for the gross wing alone is known (see Refs. 5 and 7). 


6. Numerical Example 
6.1. GENERAL 


To conclude the discussion of pressure distributions, we now consider the 
application of the preceding results to a particular example which may be of some 
practical interest. This consists of a wing of delta plan form, with rhombic cross 
sections, having a parabolic arc root section, combined with a fuselage which is also 
of parabolic section but may have a blunt base. 


(x)= 27x (1 —x) (68) 
and R(x)=2e (1 *) (wih a>l), . . . (9) 

a a 


so that the wing is of aspect ratio 4s and root thickness/chord ratio +, while the 
fuselage is of fineness ratio « and has a blunt base (if a> 1) of radius 
Ry=(2¢/a*) (a—1). 
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The fuselage area distribution is 


Sy (x)= x? (a—x)’, . (70) 
and its volume is 1 Sar 10a" — 1S$a+ 6). 


The exposed wing area distribution is 


2¢ 
where R,= sa? (a—x). 


It is convenient at this stage to introduce the parameter 


Sa 


2A represents the ratio of the fuselage slope at the apex x =0 to that of the wing plan 
form (thus when A = 4 the wing and fuselage merge at the apex). 


Then R,=24 (1 =). 
a 
so that (1- 2a) (—1+2A+ x+ 
47s a 
( 1+2A+ = x . (74) 
a a a’ 
and, from equation (17), 
Fy (x) 4 4A 72X r 440 
-6 (1-2) x+ ( - +204 4) x: 
| 12a 8A? 
The wing volume is Vy 5(1—2A)?+ (1—2A)+ (76) 
a a 


6.2. JUNCTION PRESSURE DISTRIBUTION 
Equation (44) now becomes 


ie 


(1 - (1-7 *) tog ( + + Fut (77) 
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where Gt (x)=4x(1-R,) ~R,+ (78) 
” Sw” (x) 
Sw*” (x) 
and Fy* (x)= Fw (x) 
4rs 


are given by equations (74) and (75). 
The contribution due to the fuselage is 
4:° 6x |B: x\| 16x 22x*] 


this can be written for comparison with equation (77) in the form 


Cor J 2? 4 Ox 6x? | Q x | 16x 22x? 
- 2{1- ABs {1 - 79 
reat (1 log (1 (/9a) 


4rs a a}| a a’ 


For a given value of a (that is for a fuselage with a given ratio of base area to 
maximum area), the value of C,w/(7s) depends only on the parameter A, which is 
proportional to the ratio of the maximum fuselage radius to the wing span, and on 
the wing slenderness parameter fs. The fuselage contribution C,» depends similarly 
on the fineness ratio < and on #3; the ratio of the fuselage contribution to the wing 
(plus interference) contribution is thus proportional to the ratio <*/(zs), or to the 
ratio of the fuselage volume to the wing volume. 


To provide a practical illustration of the foregoing example, the following cases 
have been computed and the results are shown in Figs. 2 to 7:- 


(a) A=4: a=1-0 (Fig. 2), | 


(Fig. 3) and 1-4 (Fig. 4) 
(b) A=}: a=1-0 (Fig. 5), 1:2 (Fig. 6) and 1:4 (Fig. 7). 
In each case the relative size of the wing and fuselage has been chosen so that 


:*/(rsa*)=2A, which is equivalent to +/s=4A. Thus typical values of +, s and « 
would be 


(a) 7=0-05, s=0-2, ¢/a=0°1 
(6) A=}: 7=0°05, s=0-4, ¢/a=0'1. 


Thecorresponding values of the ratio(fuselage volume) / (exposed wing volume)are: 


a=1-0 1-2 1-4 
3-142 7-540 14-450 
Led 0-967 1-839 2-861 
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The fuselage area distribution is 


and its volume is (10a? — 15a + 6). 


The exposed wing area distribution is 


Sw (x)=4rsx? (1 -—x) (1 -R,”). » CED 
2: 
where R,= (a—X). 
Sx sa 


It is convenient at this stage to introduce the parameter 


sa 


2A represents the ratio of the fuselage slope at the apex x =0 to that of the wing plan 
form (thus when A = 4 the wing and fuselage merge at the apex). 


Then R, 2 (1 =), 
a 
so that 2 (1 — 2A)? + 6 (1 — 2A) 1+2A+ X+ 
4rs a 
a a a 
and, from equation (17), 
=6(1—2d)( )x+ 2 5 
TS 2) 2 
The wing volume is Vy = - Isc — 2X)? + = (1 —2A)+ =| (76) 
a a’ 


6.2. JUNCTION PRESSURE DISTRIBUTION 
Equation (44) now becomes 


Cals 6x), f 4R, | 
5,” (x) log (485 (1 -8a 11 


—8A?x (1 4 (1 - lo (x) + (x)+G* (77) 


a 
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where G* (x)=4x(1-R,) (78) 
Sw” (x) 
*” (x) = 
(x) 
and Fw* (x) Fw (x) 


are given by equations (74) and (75). 


The contribution due to the fuselage is 


a}\ a a 


4:? 6x  6x° Bz x 
this can be written for comparison with equation (77) in the form 
(79a) 


E (1 6x? 
a 


16x 22x* 
4rs rsa? a 


a} a a’ 


log 3s (1 - 


For a given value of a (that is for a fuselage with a given ratio of base area to 
maximum area), the value of C,w/(>s) depends only on the parameter A, which is 
proportional to the ratio of the maximum fuselage radius to the wing span, and on 
the wing slenderness parameter 8s. The fuselage contribution C,» depends similarly 
on the fineness ratio < and on (3; the ratio of the fuselage contribution to the wing 
(plus interference) contribution is thus proportional to the ratio <*/(7s), or to the 
ratio of the fuselage volume to the wing volume. 


To provide a practical illustration of the foregoing example, the following cases 
have been computed and the results are shown in Figs. 2 to 7:- 


(a) A=4: a=1-0 (Fig. 2), 1-2 (Fig. 3) and 1:4 (Fig. 4) 
(6) A=}: a=1-0 (Fig. 5), 1-2 (Fig. 6) and 1-4 (Fig. 7). 
In each case the relative size of the wing and fuselage has been chosen so that 


2*/(7sa*)=2A, which is equivalent to 7/s=4A. Thus typical values of +, s and « 
would be 


(a) A=4: 7=0:05, s=0-2, ¢/a=0'1 
(b) A=}: 7=0-05, s=0-4, ¢/a=0'1. 


Thecorresponding values of the ratio(fuselage volume) / (exposed wing volume)are: 


a=10 1:2 1-4 
A=4 3-142 14°450 
0-967 2-861 
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(g) FUSELAGE ALONE 


Bs =0- 


TYPICAL WING 


S$=0-2 A=0'8 
WING tT =0°05 
FUSELAGE €= 0-10 


=-2r ah nc INTE RENCE 
ats (b) WING + INTERFERENCE 


Ficure 2. Wing-fuselage junction pressure distributions: 4=4, a=1-0 


—— —— —— Gross wing alone (centre section) 
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(q) FUSELAGE ALONE 


+0 


WING + INTERFERENCE P 


Ficure 3. Wing-fuselage junction pressure distributions: A=4, a= 1-2. 


—— —— —— Gross wing alone (centre section) 
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FicureE 4. Wing-fuselage junction pressure distributions: \=4, a=1°4. 
—— —— —— Gross wing alone (centre section) 
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Ficure 5. Wing-fuselage junction pressure distributions: 4=4, a= 1:0. 


— — Gross wing alone (centre section) 
Approximate method (equation (81)) 
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FiGuRE 6. Wing-fuselage junction pressure distributions: 4=4, a=1-2. 


— — — Gross wing alone (centre section) 
Approximate method (equation (81)) 
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Ficure 7. Wing-fuselage junction pressure distributions: A=4, a= 1-4. 
—— ——- —— Gross wing alone (centre section) 
- Approximate method (equation (81)) 
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The ratio of the base area of the fuselage to its maximum cross-sectional area 
takes the values 0, 0-309 and 0:666 respectively. The three examples with A=} 
thus represent combinations with relatively large fuselages and small wings, while 
those with A=} have fuselage and wings more nearly equal in volume. 


In spite of the differences in relative fuselage and wing volumes, the general 
features of the pressure distributions are very similar for the two cases A=4 and 
A=}. The parabolic fuselages are, of course, familiar and were originally investiga- 
ted by Lighthill. The pressure distributions for the fuselages alone show two 
principal features: the suction peak aft of the position of maximum radius, which 
moves rearwards and decreases in magnitude as the Mach number increases, and the 
rapid recompression following this suction peak, leading for the pointed bodies 
(a=1) to a logarithmic singularity in the pressure at the rear end.t The additional 
pressure due to the wing and interference is small compared with that due to the 
fuselage over the forward half of the wing-fuselage junction, but towards the rear it 
shows a strong negative gradient which helps to counteract the compression due to 
the fuselage. For the combinations with pointed fuselages (a= 1-0), there is still a 
resultant compressive gradient over the rear 25 per cent, which is more severe for 
the larger fuselage (A=4). For the combinations having the largest base area 
(a= 1-4), on the other hand, the adverse gradients produced by the fuselage are small 
and the overall pressure shows a strong favourable (negative) gradient up to the 
trailing edge. The combinations with a=1:2 are naturally intermediate between 
these two extremes. 


The pressure distributions, according to slender-wing theory, along the centre 
sections of the gross wings alone are also shown for comparison; it is seen that the 
addition of a fuselage causes a considerable modification in all cases. 


6.3. PRESSURE ALONG THE TOP OF THE FUSELAGE 


Equation (55) of Section 4 has been evaluated for the two cases with a sharp- 
ended fuselage (a=1,4=4 and }). The results are shown in Fig. 8. It will be seen 
that the pressure variation round the fuselage is everywhere quite small compared 
with the general pressure level, at least in this particular case. 


6.4. THE APPROXIMATE METHOD OF SECTION 5 


The method has been applied to the present example for the case \=}, which 
represents the smaller of the two types of fuselage considered. For this wing the 
root pressure distribution is” 


Cow” 


+Lighthill * points out that this singularity is simply a defect of slender-body theory and would 
not occur in a more complete potential-flow theory or, of course, in practice; but some recom- 


pression of this type does occur for pointed bodies, becoming less severe (as indicated by 
the theory) for bodies with blunt ends. 
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(b) A= 


AC, 


Ficure 8. Change of pressure round the fuselage 
AC, =C, (top) — C, (junction). 
-— — Approximate theory. 


From equation (64) it is found that, at the junction, 


ea if 3x (1+ ') log (1 - *)| (1+ =) + 


Ars a 
(81) 
a (a-X) 


The results are plotted in Figs. 5, 6 and 7 for comparison with the more accurate 
calculations of Section 6.2. It will be seen that the agreement is good over most of 
the fuselage length and that even near the apex, where R, increases to 4, the general 
shape of the overall pressure distribution is well predicted. 

A similar comparison with regard to the change of pressure round the fuselage 
is made in Fig. 8. In this case (A=}, a=1-0) z, is a constant multiple of R so that 
the last term in equation (66) vanishes and we have simply 


AC, = + Sar” log 2. 


The agreement with the more accurate equation (55) is again quite satisfactory, 
although, since the value of AC, given by (55) is itself everywhere small, the shape of 
the curve is not very well predicted by the approximate theory; but the actual 
numerical difference between the two results never exceeds 7s, and this is as good 
as for the pressure at the wing-fuselage junction. 


May 196] 125 


irea 
hile 

ac, 
eral 
and 
iga- 
Iwo 
ich 
lies 
nal 
the =| | 
r it 
to 
la 
for 
ea 
all 
he 
en 
re 
he 
n 
d 
h 
e 

3 
a 


C LOCK 


1-2 asl-O 


aal-O 
NEWBY WING 
xe) 


0-6 


0-4 asl-4 


0-2 
WiING-FUSELAGE COMBINATIONS 

FUSELAGES ALONE — 


| | 


Ficure 9. Wave drag factors K,: A=. 


6.5. OVERALL WAVE DRAG 


The appropriate extension of the Ward-Lighthill slender-body theory” to wing- 
fuselage combinations of the general type considered in this paper is given in Ref. 4. 


This method has been applied to the present example, and the results, in the form of 
the factor (cf. Ref. 8) 


= 
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Figure 10. Wave drag factors K,: 


are shown plotted against @s in Figs. 9 (A=4) and 10 (A=}). For comparison we 
have also included values of K, for the three fuselages (a=1-0, 1-2 and 1-4) alone, 
and for some typical wings alone: the gross wing used in this example (“Newby” 
wing”) and Lord’s® wings 4 and 5. It will be seen that in both cases (A=4 and }) 
the combinations with fuselages pointed at both ends (a=1-0) have drag factors 
considerably lower than for the gross wing alone; for the combination with the 
larger fuselage (A= 4) the values are close to Lord’s wing 5 while with the smaller 
fuselage (A = }) they are close to Lord’s wing 4. These low values are largely due to 


May 1961 127 


0-8 
me 


the favourable interference produced by the compression near the rear of the 
fuselage acting on the backward facing surface of the wing; since this compression 
will not be so large in practice, the drag results may be over-optimistic. 


In both these cases the drag factors are also lower than for the corresponding 
fuselage alone over a wide range of values of fs, but this is no longer true for 
combinations with blunt-based fuselages (a= 1-2 and 1:4); the combinations with the 
larger fuselages (A= 4) have values of K, very similar to the fuselages alone, while 
with the smaller fuselages (A = }) the values of K, are considerably higher. Blunting 
the base of the fuselage causes in all cases a considerable reduction in K,; whether 
any practical use can be made of this depends on the extent to which such a base 
could be incorporated into the engine design. 


7. Conclusions 


The method described in the present paper provides a means of calculating the 
pressure distribution and wave drag in supersonic flow for a wide class of 
symmetrical wing-fuselage combinations at zero lift. For the method to be strictly 
applicable, the wings should be “thin” in transverse cross section, both in themselves 
and relative to the fuselage. It is now known that in practice the first of these 
conditions can often be relaxed; it remains to be seen how far this is possible with 
regard to the second. If the wing is not thin relative to the fuselage, it might be 
desirable to replace the equations 


Sw (x)=4 Zw (x, dy 


R 
Swr (x)=4 [ Zw (x, 9) dy 
by the appropriate exact expressions, taking into account the curvature of the 


fuselage inside the wing, although it is unlikely that this would affect the accuracy 
significantly. 


A further generalisation of the method would be to combinations in which the 
fuselage extends ahead of the wing root leading edge; these would present no 
additional difficulties provided that the overall area distribution remained sufficiently 
smooth. It might also be useful in estimating sting interference for all-wing models 


in which a sting of circular cross section is blended into the wing shape near the 
trailing edge. 


The principal results of the paper can be summarised as follows. 
7.1. OVERALL PRESSURE DISTRIBUTION 


The pressure coefficient is given by equations (11) to (15) in the form C,p+ C,w. 
where C,» is the contribution of the fuselage alone (equation (12)) and C,w is the 
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contribution of the wing and interference effects (equations (13) and (15)). For a 
general position on the wing there are two integrals in equation (15) which have not 
been evaluated explicitly; it is suggested that this could best be done by purely 
numerical methods. 


7.2. WING-FUSELAGE JUNCTION PRESSURES 


Along the wing-fuselage junction, C,w is given by equations (36) or (36a), which 
are valid for wings of quite general shape; if the wing thickness distribution is known 
as a polynomial & a, (x) y" in », the function G (x)= G, (x) can be obtained from 
equation (37). Explicit formulae for C,w!; for wings of rhombic or parabolic cross 
section are given by equations (44) and (50) respectively. 


When the ratio R, of the fuselage radius to the local semi-span is fairly small, 
a considerable simplification may be made and the following formula results : 


where C,,» is the pressure on the fuselage alone (equation (35)), C,w°” is the pressure 
along the centre line on the gross wing alone (equation (62)), and C,, is the 
interference pressure given by 


- log (=) +Fwe—4R . (64a) 


d {(z,RY)\ 
dx | R | 


this is valid for a wing of general thickness distribution, and will be particularly 
valuable if the pressure C,w°”’ due to the gross wing alone is already known. 


7.3. PRESSURE ALONG THE TOP OF THE FUSELAGE 


The change of pressure round the fuselage, from the wing junction to the top, is 
given by equation (54) in general; for wings with rhombic or parabolic cross sections 
this takes the form (55) or (56) respectively. The corresponding approximation when 
R, is small is given by equation (66). In the examples considered in Section 6, this 
change of pressure was relatively small in all cases. 


7.4. EXAMPLE 


A simple example has been considered in Section 6. consisting of a “Newby” 
wing” combined with a parabolic body of revolution (which may have a blunt base). 
The drag results (Figs. 9 and 10) show favourable interference between the wing and 
the fuselage: this is due to the compression near the wing trailing edge produced by 
the fuselage (see Figs. 2-7). When the fuselage is pointed at the rear the theoretical 
pressure distributions have a spurious singularity there which may invalidate the 
drag results, but when the fuselage is blunt the pressure gradients are favourable 
over almost the whole of the wing and so could be realised in practice. 
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The Axial Decay and Radial Spread of a 
Supersonic Jet Exhausting into Air at Rest 


J. H. FRAUENBERGER* and J. G. FORBISTERt 


SuMMARY: The distribution of total pressure, stagnation temperature and gas 
velocity is determined for the subsonic region of a supersonic jet emerging from 
a solid propellant rocket motor into air at rest. Measurements of Pitot pressures 
in the supersonic flow region have shown that the flow follows the axis of an 
inclined nozzle within the order of accuracy of the measurements. Nozzle cone 
angle, and probably the expansion ratio, affect the jet distribution significantly. 


1. Introduction 


Knowledge of the effect of the jet from a rocket motor on remote surfaces such 
as the control surfaces of a missile, the launcher, the launching aircraft or other 
structure is often desirable for design purposes. Some firings of a standard solid 
propellant rocket motor were arranged to provide information for estimating effects 
of this kind. Total pressures and total temperatures were measured at a number 
of positions in the jet. Velocities of the gas-air mixture have been computed from 
these measurements. 


Investigations on free subsonic jets (length of core, quantity of entrained air, 
and decay of velocity) have been made by Kuethe’?’ and Wuest®. The present 
investigation into the axial decay and radial spread of a supersonic jet exhausting 
into still air considers both the early supersonic régime and the later subsonic jet. 


NOTATION 
a_ velocity of sound 
ad critical velocity of sound 
c gas velocity (ft./sec.) 
diameter (in.) 
y=c,/c,, ratio of specific heats 


M*=c/a’, Mach number 


*Now in Germany 
+Atomic Weapons Research Establishment. 
Aldermaston 


Formerly Rocket Propulsion Department, 
Westcott. 


—- 
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M=c/a, Mach number 
m mass flow rate 
P pressure (Ib./in.? absolute) 
pressure (Ib./in.*? gauge) 
P* Pitot pressure (Ib. /in.? absolute) 
p* =(P* — P,), Pitot pressure (Ib./in.* gauge) 
R_ gas constant 
r radius (in. or ft.) 
T temperature (°K) 


x distance of point of consideration from nozzle exit (in. or ft.) 


Suffixes 
G gas 
A air 
c¢ combustion chamber 
0 stagnation value 
m gas-air mixture 
e nozzle exit 


t nozzle throat 


2. The Development of the Jet 


Before describing and discussing the tests and their results it may be of some 
interest to consider the development of the gas stream and the use of Pitot tubes for 
measurements of total pressures in supersonic and subsonic flows from a purely 
theoretical aspect and to give a general picture of the results. 


2.1. PRESSURE MEASUREMENT BY A PITOT TUBE 


On placing a Pitot tube into a jet stream of supersonic velocity c, a compression 
shock will form ahead of it and a sudden change in pressure, velocity and density of 
the gas will occur. The measured pressure P* (pressure after the shock at the 
stagnation point) is less than the reservoir or total pressure P, corresponding to the 
velocity c. In a subsonic stream the Pitot tube measures the total pressure. 


When the total pressure, P,, of the supersonic flow is known, the Mach number 
is given by the formula 


pr y+1 
M*? y+ 1 
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In this formula M*=c/a’ is the Mach number related to the critical velocity of 
sound a’, and y is the ratio of the specific heats of the gas. This relation is convenient 
for the early supersonic flow from a nozzle before shocks develop. where the total 
pressure is equal to the chamber pressure of the rocket motor, that is P,=P.. Else- 
where in the flow knowledge of the free stream static pressure, P, is required to 
determine the Mach number, which is then given by 


aes 
. | forM>1 
y+1 y+1 


where M =c/a is the Mach number related to the local velocity of sound. 


The dependence of M* on P*/P, and of M on P*/P for y=1:25 and y~-14 
are shown in Fig. | and Fig. 2 respectively. 

No other measurement than P* has been made in the work described in this 
report, which is mainly concerned with the flow some distance from the nozzle. It is 
assumed that the static pressure in the subsonic portion of the jet is equal to the 
atmospheric pressure. 


2.2. THE INFLUENCE OF AN UNDER-EXPANDED NOZZLE ON THE PRESSURI 
DISTRIBUTION IN THE JET STREAM 


For a parallel jet leaving an under-expanded two-dimensional nozzle and 
expanding into still air, Prandtl gives a wave diagram which is reproduced in 
Fig. 3. In this simple case the jet, as it flows from the nozzle, expands and contracts 
between the pressures P, and P, where 


P, is the initial stream pressure and P, is ambient pressure. In practice, this simple 
cyclic flow would gradually break down owing to the development of shock waves 
and mixture with the surrounding air. 


FicureE 3. Wave diagram of supersonic jet from under-expanded parallel nozzle. 
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Ficure 4. Jets from under-expanded (a) and fully expanded (4) nozzles 


The flow from an under-expanded axially symmetrical parallel or conical nozzle 
will also have alternate regions of high and low pressure with, however, a more 
complex pattern of shock waves. A large number of detailed photographs of the 
air flow from conical nozzles with various expansion ratios and cone angles have 
been taken by Fraser and Rowe.’ Photographs taken in the course of the present 
work are shown in Figs. 4 and 5. Fig. 4(a) shows the jet from the under-expanded 
nozzle and the periodic shock zones of decreasing intensity are clearly visible. 
Similarly, Fig. 4(b) shows the weaker shock zones in the jet from a fully expanded 
conical nozzle. In this case the shocks are not so strong, but nevertheless do exist, 
because the flow from the conical nozzle is not parallel and because of viscosity 
effects. 


2.3. GENERAL STRUCTURE OF THE JET 


A jet leaving a nozzle at initially supersonic velocity contains along its axis the 
gradually narrowing core of the original gas stream surrounded by an annular region 
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(a) 38°75 in. from under-expanded nozzle. 


(b) 23-5 in. from fully expanded nozzle. 


FicureE 5. Pressure measurement in jets from under-expanded and fully expanded nozzles. 


where gas and entrained air have intermixed. As mixing extends across the cross 
section, temperature and velocity fall and further mixture of air with gas and already 
entrained air takes place at subsonic speed (Fig. 6). 
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FiGcure 6. General structure of a jet 


It will be assumed for the subsonic portion of the jet that the static pressure 
equals the ambient pressure. Since, at the beginning of the subsonic portion, M 
equals unity, the supersonic Pitot tube formula will reduce to 
px 
p 


For y=1:25 (which is an average value for the products of combustion of our 
motor), P*=1-8P and for air, P*=1:893P. Assuming a static pressure of one 
atmosphere and that the product gases are air rich by the time they reach the sub- 
sonic region, the subsonic portion begins with a Pitot pressure of 27:5 Ib./in.* 
absolute. 


If the gases leaving the nozzle are fuel rich, afterburning is likely to occur, 
which may have some influence on the physical properties of the mixture. 


3. Test Apparatus 


The solid propellant motor used in the present experiments burned for some 
2:5 seconds at a chamber pressure p. of 1,000 to 1,100 Ib./in.*, as can be seen in 
Fig. 7. The mean mass flow rate was 33 Ib./sec. Assuming a chamber pressure of 
1,030 Ib./in.* absolute and a gas temperature of 2,473°K, the exit velocity for 
isentropic expansion to atmospheric pressure is 7,380 ft./sec. at a gas temperature 
of 1,055°K. The gas composition is shown in Table I. 


Owing to the arrangement of the test bed, and to reduce interference between the 
jet and the ground, inclined nozzles were used on the motor. Nozzles having inclina- 
tions of 15°. 11° and 7:5° were selected from a range of existing types. The 7-5° 
nozzle restricted the range for the measurements to a distance of some 40 ft. from 


TABLE I 
GAS COMPOSITION IN MOTOR 
Component H,O CO, CO H, OH H N, NO 
Percentage by volume 20 8:96 43-39 15-23 0-30 0-12 12-26 0-01 
Pressure — 70 atmospheres Temperature = 2473 K 
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FIGURE 7. Pressure-time record for motor. 


the exit plane for reasonably symmetrical entrainment of air, assuming a cone angle 
for the jet of 20° (Fig. 8). The 15° and 11° nozzles were used for measurements in 
the supersonic and the earlier part of the subsonic regions respectively. 


In most of the experiments nozzles having an included angle of 30° were used, 
but some tests were made with nozzles of 15° included angle for comparative 
purposes. All had a throat diameter of 2°55 in. and most had a ratio of exit area to 
throat area of 4:1, i.e. they were under-expanded, giving an exit pressure of 
41 lb./in.* absolute. Some nozzles with an expansion ratio of 7:7: 1, however, were 
used; these were almost fully expanded, with an exit pressure of 16:3 Ib./in.* 
absolute. The corresponding Mach numbers M at the nozzle exits were 2°69 and 
3-22 respectively. 


Pressure measurements in the jet were made by Pitot tubes and two types proved 
reliable and able to withstand the gas temperatures after some development. One 
was the water-cooled tube shown in Fig. 9(a), consisting of the Pitot tube itself within 
a larger supporting tube which allowed water to flow down the annular space for 
cooling purposes. The other consisted of a plate (Fig. 9(b)), drilled with 5 or 16 
holes according to the size of the jet. These were employed in the supersonic region. 
In the sonic and subsonic region single tubes, as described but uncooled, were 
secured at 8 in. intervals on the vertical portion of an instrument carriage (Fig. 9(c)) 
that could be inclined at different angles in accordance with the inclination of the 
motor nozzle. The single tubes were provided with hemispherical heads of diameter 
six times the bore of the internal tube; likewise the plate type had a rounded front 
edge, the plate thickness being six times the bore of the holes. No attempt was made 
to move these devices during a firing, either radially across or axially along the jet, 
owing to the short duration of each firing and the slightly changing chamber 
pressure. 
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Ficure 8. Spread of jet from a 7°5° inclined nozzle assuming a jet angle of 20°. 
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Ficure 9. Measuring arrangements. 
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The Pitot tubes were aligned optically with the motor nozzle. A hollow Tufnol 
body, which contained a weak light source giving a beam through two fine orifices, 
was machined to fit in the nozzle; by rotating this body the errors of manufacture 
were eliminated and the jet axis could be located in relation to the Pitot tubes to 
within }°. The alignment was liable to changes during the run in accordance 
with the extent of erosion of the nozzle. 


The pressures were recorded electronically (Miller recorder) by means of a 
strain-gauge pressure transducer, consisting of a short elliptical tube with strain 
gauges bonded to each end of the major and minor axes. On pressurisation of the 
tube the resistance of one pair of gauges increases and that of the other decreases; 
the gauges were connected to a Wheatstone bridge to give the maximum sensitivity. 


Temperatures were measured by chromel-alumel and platinum-platinum/ 
rhodium thermocouples fixed between the Pitot tubes (Fig. 9(c)). 


In addition to the instrumentation, cinéfilms of the rocket jet were taken by a 
high speed camera remotely operated at 24 or 48 frames per second. 


4. Measurements in the Region of Supersonic Flow 
4.1. GENERAL 


The values of the Pitot pressures were taken from the records at times 
corresponding to a pressure of 1,030 Ib./in.? absolute (1,015 Ib./in.* gauge) in the 
chamber, because conditions were steady at this pressure for about 4 second. All 
these pressures are plotted in Figs. 10, 11 and 12 and their values on the axis are 
given in Table II. In general the results confirm the cyclic nature of the flow in the 
supersonic region. The variation of Pitot pressure across a section of the jet has 
two forms which alternate along the jet. In one the Pitot pressure is high at the axis 
and falls rapidly towards the boundary of the jet. This corresponds to a high static 
pressure at the axis falling to ambient pressure at the boundary. In the other case 
the Pitot pressure is low at the axis, rises to a maximum near the boundary and falls 


TABLE II 


MEASURED PITOT PRESSURES ALONG THE AXES OF JETS EMERGING FROM UNDER-EXPANDED AND 
ALMOST FULLY EXPANDED NOZZLES 


Both Nozzles: included angle= 30°, combustion chamber pressure, P,= 1030 Ib./in.? absolute 
Under-Expanded Nozzle 
area ratio=4:1, calculated Pitot pressure at nozzle exit, P* =373 lb./in.2 absolute 
Distance from nozzle x (in.) 28-5 38°75 104-5 109-5 157 
Pressure p* (Ib./in.* gauge) 180 75 23 


Fully Expanded Nozzle 
area ratio=7'7:1, calculated Pitot pressure at nozzle exit, P*=210 Ib./in.2 absolute 


84:5 


Distance from nozzle x (in.) 12 23-5 4 32-5 


Pressure p* (ib./in.* gauge) | 30 297 330 120 125 
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Ficure 10. Pitot pressure distribution across the jet from under-expanded nozzle. 
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Figure 11. Pitot pressure distribution across the jet from fully expanded nozzle. 


sharply to ambient pressure. This corresponds to a low static pressure at the axis 
rising to atmospheric pressure at the boundary. The levels of the Pitot pressure fall 
from cycle to cycle due to losses caused by shocks and mixing with air. 


Referring back to Fig. 4(a), the jet from an under-expanded nozzle with a cone 
angle of 30° may be seen. The expansion of the flow from the nozzle exit is clearly 
shown and the shock zones correspond to regions of high pressure on the axis of the 
jet. (Compare Fig. 4(b) for the jet from a correctly expanded nozzle.) 
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4.2. EFFECT OF THE NozzLE Exit CONE ANGLE 


Measurements of the pressure distribution across the rocket jet were made at a 
distance of 104 in., or 41 throat diameters, from the nozzle exit, using under- 
expanded nozzles with exit cone angles of 30° and 15°. The nozzles were inclined 
upwards at 15°. 


The results are plotted in Fig. 12, where position is measured from the axis of 
the nozzle. They show that the jet follows the axis of the inclined nozzle to the 
accuracy of the experiments, say within +}° (+4 in. at 104 in.). To this accuracy 
the cone angle of the nozzle does not affect the deflection of the jet. 


ANG 
80 
| 
p* (LB./IN.? | | 
GAUGE) | : 
| NN 30 INCLUDE NG 
| 
\ | 
| 
u x 
20 ASH 
2 | 
4 | 
3) + 4 5 4 10 2 


FiGureE 12. Pressure distribution across the jet from under-expanded nozzle 
with different included angles. Pressure at 104 in. from nozzle. 


The smaller nozzle cone angle does, however, result in a more concentrated jet 
at this distance. The Pitot pressure (gauge) falls to half its peak value on the axis at 
radial distances of 2:5 in. and 3-5 in. for cone angles of 15° and 30° respectively. 
Radial distances from the axis may be expressed more conveniently as angular 
deflections from the axis referred to a centre at the exit of the nozzle. Then the 
angular deflections, for which the Pitot pressure is one half of its peak value on the 
axis at 104 in. distance, are 14° and 2° for the nozzles with cone angles of 15° and 
30° respectively. 


4.3. THE FLOW NEAR THE NOZZLE AND THE EFFECT OF NOZZLE EXPANSION RATIO 


For comparative purposes extensions were welded to the expansion cones of 
some nozzles to give an area ratio of 7:7, which corresponds to about full expansion 
for a chamber pressure of 1,000 Ib./in.* The effect of area expansion ratio will be 
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FiGure 13. Pressure measurement 12 in. from fully expanded nozzle. 


greatest in the initial portion of the jet, so measurements of Pitot pressure were made 
in the early supersonic portion of jets from nozzles having a cone angle of 30° and 
area ratios of both 4 and 7.7. 


The results for the flow from the under-expanded and fully expanded nozzles 
are given in Figs. 10 and 11 respectively. The letters (a) to (f) refer to different 
distances from the nozzle exit; they are also shown in Fig. 6 which presents the 
general structure of a jet emerging from an under-expanded nozzle. It should be 
noted that the scale in Fig. 6 is approximately true only for the first and second wave 
of the jet. The pressure distribution 12 in. from the exit of the fully expanded nozzle 
(Fig. 13) is shown in Fig. 11. There are no shocks in the flow before this point so 
that the total temperature and pressure are known. The Pitot pressure, therefore, 
determines the flow at this section completely. The Pitot pressure of about 30 Ib. /in.* 
gauge on the axis at this distance indicates a static pressure of about 2 Ib./in.? 
absolute and this is very close to the value 3 Ib./in.* absolute corresponding to 
conical flow from the nozzle. Based on the static pressure of 2 Ib./in.*, the gas 
velocity is 8,235 ft./sec. (M=4°5). The flow from the under-expanded nozzle is 
similar (Fig. 10). The steep and variable pressure gradient across the flow did not 
permit of accurate determination of the spread of the jet in this region of the flow. 
The spread from the under-expanded nozzle is, however, comparatively slow, 
amounting to about + 12 in. at 104 in. from the nozzle. This is +6°5’. 


5. Measurements in the Region of Subsonic Flow 


5.1. GENERAL 


All tests for measuring pressures and temperatures in the subsonic region were 
made with under-expanded nozzles of 30° included angle. The beginning of the sub- 
sonic portion was estimated from the results to be at a distance of 17 ft. from the 
nozzle, All Pitot pressures were taken from the records at a time when the chamber 
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pressure in the motor was 1,030 Ib./in.* absolute. As previously stated, cinéfilms at 
24 and 48 frames per second were taken of many of the firings. In all cases after- 
burning occurred intermittently and appeared as a short flame within the range of 
13 to 20 ft. from the nozzle. The flame was recorded on a single frame about once 
every five frames (Fig. 14). The thermocouples did not respond to oscillations of 


Ficure 14. Intermittent afterburning during a firing. 


this frequency and the temperature record was smooth, indicating a value between 
the temperature with afterburning and that without afterburning. The recorded 
temperature would consequently be heavily weighted towards the value for no 
afterburning by a factor of at least five. 


Assuming no afterburning and no heat losses by conduction and radiation, we 
have 
MeCyaT + = (Me + My) ComT m- 


where c, is the specific heat at constant pressure and the suffixes are defined in the 
Notation. Hence it is p~ sible to correlate a certain mixture ratio and value of y to 
a measured stagnation temperature; the gas velocity can be computed using 


P* 
5.2. AXIAL DECAY 


The pressures, temperatures and velocities along the axis of the jet are shown 
in Fig. 15. In the non-dimensional scales used for the same gas properties in 
Fig. 16 any reference to the conditions at the under-expanded nozzle was avoided 
by taking the throat diameter d,, throat velocity c,, motor pressure P. and motor 
temperature 7,; the average air temperature was + 10°C. 


5.3. RADIAL DECAY 


At each distance from the nozzle indicated in Fig. 15, Pitot pressures and 
temperatures were measured at a number of distances from the axis of the subsonic 
jet by Pitot tubes and thermocouples mounted on a stand, illustrated in Fig. %c). 
Individual observations have not been reproduced but have been analysed assuming 
that the radial decay in the subsonic jet is a function of the angle from the axis of 
the nozzle only and is independent of the distance from the nozzle. 


The angular displacements at which the pressure, temperature and velocity 
fall to 3/4, 1/2, 1/4 and 1/10 of their values on the axis are given in Table III 
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on p. 146. Complete diagrams of the axial and radial decay of temperature, 
pressure and velocity are given in Fig. 17. 


5.4. COMPARISON WITH ANOTHER INVESTIGATION 


While this paper was being prepared a report became available of an American 
investigation in which an experimental study was conducted of the spreading and 
axial decay of free supersonic jets exhausting into quiescent air. Supersonic heated- 
air jets and three types of solid propellant rockets were used, the Mach number M 
varying between 1-4 and 3:53. In the following comparison only the three solid 
motors will be considered. 


The nozzles were of the convergent-divergent type with included angles of the 
exit cone of 30° to 60°, exit diameters of 2:0 to 2:59 in. and Mach numbers of 2°82 
to 3:53. Although not actually stated, it will be assumed that the nozzles were 
fully expanded and that the ratio of specific heats was 1-25. The maximum distance 
at which measurements were made was 26 ft. approximately. 
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Ficure 15. Pitot pressure p* (lb./in.? gauge), stagnation temperature T,, (°C) 
and velocity c (ft./sec.) along the axis. 
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Ficure 16. Pitot pressure, stagnation temperature and velocity along the axis 
expressed in non-dimensional ratios, 
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FIGURE 7. 
Top diagram: Distribution of stagnation temperature in jet stream. 
Middle diagram: Distribution of Pitot pressure in jet stream. 
Lower diagram: Distribution of velocity in jet stream. 


The results of Ref. 5 indicate lengths of the supersonic jets (varying between 
6°25 and 9 ft.) which are much smaller than that from the motor used for the present 
work (for which it is 17 ft.) and that the axial pressure and temperature decay starts 
much earlier; e.g. at a distance of x=26-°5 ft. our measured pressure is about 12 
times and our temperature is about 5 times larger than the respective values of 
rocket A (M = 3-53) in Figs. 18 and 19. 


TABLE III 


THE ANGULAR DISPLACEMENT FROM THE JET AXIS FOR GIVEN DECAY RATIOS OF TEMPERATURE, PRESSURE 
AND VELOCITY IN THE REGION OF SUBSONIC FLOW FROM AN UNDER-EXPANDED NOZZLE 


Temperature Pressure Velocity 

Ratio Angular Displacements 

3/4 1-8 1-7 

1/2 3-5 3°2 4 

1/4 3°2 4-7 6 

1/10 6°4 
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Ficure 19. Stagnation temperature along the axis compared with values from Ref. 5. 
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A more adequate comparison of the axial pressure and temperature decay in 
the subsonic jet can be made by introducing non-dimensional parameters. Since 
the throat diameters of the rockets mentioned are not known, the values found in 
the present work have been recalculated in terms of the ratio of distance from the 
nozzle to nozzle exit diameter x/d,. The results are plotted as curves D in Figs. 18 
and 19; the pressure at the exit is 41 Ib./in.* absolute and M =2-69 at the exit. The 
results of Ref. 5 are plotted as curves A (M =3-53), B (M =3-42) and C (M =2°82). 


When comparing these curves it should be noted that curve D applies to an 
under-expanded nozzle of cone angle 30°, whereas A, B and C apply (by inference) 
to fully expanded nozzles of cone angles 50°, 60° and 30° respectively. If the 
results of curve D are replotted in terms of the theoretical exit diameter for full 
expansion to atmospheric pressure (p= 14-7 Ib./in.? absolute and M =3-28), curves 
E of Figs. 18 and 19 are obtained. It is clear that the results are not compatible. 
Further experiments are necessary to establish the effects of nozzle expansion ratio 
and cone angle. 


(Curve C has been omitted from Fig. 19 since, in contrast to the curves for 
rockets A and B, this rocket produced a jet with spontaneous afterburning, which 
shifted the temperature scale downstream.) 


6. Conclusions 


(i) Observations in the supersonic flow unaffected by disturbances from the 
exit of a conical nozzle (included angle 30°, area ratio 7°7:1) are consistent with 
very nearly conical or “source” flow at the exit plane of the nozzle. 


(ii) The flow from a nozzle of area ratio 4:1 inclined at 15° from the motor 
axis follows the axis of the nozzle within +}°. This is true for nozzles having 
included angles of 15° and 30°. 


(iii) The flow from a nozzle of area ratio 4:1 with 15° included angle is more 
concentrated than that from one with 30° included angle. At about 40 throat 
diameters from the nozzle the Pitot pressure falls off to one half its value on the 
axis at angular displacements of 14° for the 15° cone and 2° for the 30° cone. 


(iv) In the region of subsonic flow the total pressure, stagnation temperature 
and velocity decay logarithmically along the axis. The results for the motor used 
are given non-dimensionally in Fig. 16. 


(v) For the motor (with under-expanded nozzle) used, the total pressure 
(gauge), stagnation temperature (above ambient temperature) and velocity fall off 
to one half their values on the axis at angular displacements of 3:2°, 3:5° and 
4° respectively. 


(vi) From the few measurements made in the region of supersonic flow and 
from comparison with the results of Ref. 5, it is clear that the angle of the nozzle 
cone, the expansion ratio of the nozzle relative to expansion to atmospheric pressure, 
as well as the combustion chamber pressure and temperature, affect significantly 
the distribution of pressure, temperature and velocity throughout the exhaust jet. 
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Compression Tests on Thin-Walled Cylinders 


M. HOLMES, B.Sc., Ph.D., A.M.L.C.E. 


(Department of Civil Engineering, University of Bristol) 


SUMMARY: The paper describes compression tests on eight thin-walled cylinders 
of 3 ft. diameter and 0-035 in. wall thickness made of aluminium alloy plate. 
The lengths of the cylinders were either 6 ft. or 9 ft. Three of the cylinders 
were tested under axial compression up to buckling failure, and the initial 
buckling load, failing load and mode of buckling were observed. A further three 
cylinders were similarly tested, but these cylinders were subjected to internal 
pressure before applying the compressive load. The internal pressure had a 
twofold strengthening effect on these cylinders. Firstly, it induced a tensile 
pre-stress along the axis of the cylinder and, secondly, it resulted in the value of 
the compressive stress at which buckling occurred being greater than the 
buckling stress value for the unpressurised cylinders. The six cylinders were 
tested in a manner which allowed the end face of each cylinder to rotate about 
a diametral axis. A parallel platen device was used in testing the last two 
cylinders (in an unpressurised condition) which restrained rotation of the end 
faces of the cylinders. These tests enabled the effect of end restraint to be 
studied, and also enabled measurements of load-carrying capacity at large axial 
deflections to be made. The initial buckling loads, failing loads and modes of 
buckling observed in the tests were compared with existing large deflection 
theory. 


1. Introduction 


It is well known that thin tubes in compression commonly buckle in one of two 
modes; for very thin tubes, diamond pattern buckles (more or less elastic) form with 
their diagonals longitudinally and circumferentially disposed, while for rather thicker 
tubes circumferential wrinkles (plastic) arise to form circumferential corrugations. 
It is also generally known that the compression loads at which such deformations 
occur are substantially below (roughly one third of) the buckling loads predicted by 
small deflection theory. 


The production of accurate thin tube for missile purposes has made available 
thin-walled cylinders with diameter/ thickness ratios higher than hitherto practicable: 
this, with the possibility of, and current interest in, pressurisation of such tubes, 
rendered compression experiments on thin large cylinders very desirable. 


By arrangement with Mr. H. B. Howard of the Ministry of Supply (now 
Ministry of Aviation), eight 3 ft. diameter tubes, some 6 ft. long and some 9 ft. long, 
were specially constructed by Saunders-Roe Ltd. This paper describes the tests 
carried out on these cylinders in the Department of Civil Engineering of the Univer- 
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TABLE I 
Buckling Load (tons) 
Specimen _ End Length (Ib./in.?) Total Taken by | Taken by 
Condition (ft.) : di 
water cylinder 
at top at bottom pressure 
| ball 6 0 0 7-9 0 7-9 
2 ball 6 0 0 5-9 0 5-9 
3 ball 6 28-4 31-0 28-8 12-9 15-9 
4 ball 9 0 0 52 0 5-2 
§ ball 20-7 25-4 9-4 16-0 
6 ball 9 8-1 12-0 15-8 3-7 12-1 
7 platen 6 0 0 8-6 0 8-6 
& platen 8-5 0 0 6-5 0 6°5 


sity of Bristol. A large testing frame available in the laboratories there was specially 
modified for the purpose and a brief description of this is included. 


The large deflection theory of thin tubes in compression, initiated by von 
Karman and by Leggett in the early part of the Second World War, has since been 
developed in the United States to apply to pressurised as well as unpressurised thin 
cylinders. At the end of this paper the results of the tests described here are 
compared with the predictions of this theory. 


Table I shows the main results obtained during the testing of the specimens and 
the photographs of Fig. | illustrate the types of buckling failure which were 
observed. 


NOTATION 
R_ radius of cylinder 
length of cylinder 
t wall thickness 

P internal pressure 

E modulus of elasticity 

f-- buckling stress for internal pressure p 
fer’ = fer (R/ Et) 
feo buckling stress for internal pressure p =0 
foro’ = (R/ Et) 

p’=(p/E)(R/t? 


2. Description of Specimens 


The specimens were all formed from L-70 aluminium alloy in 20 s.w.g. sheet. 
The 6 ft. long specimens were manufactured with one longitudinal lap-welded seam 
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(a) Specimen 1. Failed at 7:9 tons, 
pressure zero. 


(b) Specimen 6. Failed at 15-8 tons, 
pressure 8:1 Ib./in.?. 


(c) Specimen 5. Failed at 
25:4 tons, pressure 20-7 
Ib. /in.?, 


(d) Specimen 3. Failed at 
28:8 tons, pressure 28-4 
Ib. /in.?. 
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For the 9 ft. specimens, circumferential lap welds were used to join three 3 ft. 
sections, the longitudinal seams in each section being staggered. In manufacture one 
of the 9 ft. specimens was damaged; after cutting away the damaged section the 
finished specimen was 8 ft. 6 in. long. 


Each specimen was fitted with end rings formed by bending a length of angle 
section to a suitable radius. One leg of this angle section was riveted to the cylinder 
and the other leg was bolted to a 3 ft. 6 in. x 3 ft. 6 in. x 4 in. steel end-plate. During 
testing, both the steel end-plates were backed by a 3 ft. 6 in. x3 ft. 6 in. x4 in. 
reinforced concrete slab and a layer of cement mortar was applied between the slat 
and steel end-plates to ensure overall contact. The various features described can 
be seen in the photographs of Fig. 1. 


3. Testing Frame 


The specimens were tested in a large frame in the Civil Engineering Laboratories 
at Bristol University. Essentially this testing frame is a strong structural framework 
built up from standard steel sections. A clear space 20 ft. high x 12 ft. square is 
provided in the framework for the accommodation of test specimens. At the base 
of the framework are beams which may be adjusted to support any desired test 
specimen. Load is applied to a specimen by a single 100-ton hydraulic ram and the 
framework is so designed that this ram may be placed in any desired position, either 
in elevation or plan, within the framework. 


Load measurements are not taken directly from the ram pressure, because of the 
relatively high frictional forces which exist in the ram. Instead, loads are measured 
by an Elliott (electrical resistance) Load Cell placed in series with the ram. 


4. Parallel Platen Device 


Specimens 1-6 were tested with a ball-and-socket joint between the loading ram 
and the reinforced concrete backing slab. In testing Specimens 7 and 8 it was desired 
to prevent rotation of the end faces of the specimens; consequently a parallel platen 
device was used for load application. This device consists of a tubular guide, 7 in 
diameter and 18 in. long, which is held rigidly under the loading ram by beams 
connecting into the structure of the testing frame. Through this guide passes a solid 
circular shaft, of the same diameter as the guide, with a mushroom end of 24 in. 
diameter. Load is thus applied from the loading ram through the shaft with its 
mushroom end to the concrete backing slab of the specimen. In this way rotation 
of the end face of the specimen is prevented. The parallel platen device is illustrated 
in Fig. en 


5. Application of Internal Pressure 
It has been mentioned already that some of the cylinders were tested under 
axial compression with internal hydrostatic pressure. The method by which this 


pressure was applied and measured is illustrated in Fig. 3. 
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(a) Specimen 8. Parallel platen test. 


(b) Parallel 
platen device 


FIGURE 2. 


As there was a considerable leakage of water at the welded seams and the end- 
plates, a procedure had to be adopted to keep the pressure constant during the test. 
Valve A was set to give a supply from the mains which just exceeded the losses due 
to leakage; this setting of valve A was left unaltered during the course of the test. 
During the test, water was bled from the system through the fine control valve B, so 
that the levels in the mercury manometer at C and D remained constant. In this way 
the pressure in the cylinder could be maintained within +2 per cent of the desired 
value during the entire course of a test. 
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FIGURE 3. Pressure system. 


MERCURY 


The use of water as a pressurising agent meant that there was a variation in 
pressure from the top to the bottom of a test cylinder. This was considered 
preferable to the use of air and the consequent dangers of a blow out. 


6. Test Procedure 


Before beginning the loading tests, measurements were made of the cylinder 
wall irregularities; where appropriate, these measurements were repeated after the 
application of internal pressure. 


In carrying out tests in which internal pressure was to be applied, the full 
pressure was applied first, then the axial compressive load was applied and increased 
in small increments up to failure. For some, but not all, of the test cylinders, 
measurements of the longitudinal shortening of the cylinder were made. Readings 
were taken with 1/1000 in. dial gauges across the four corners of the end-plates; 
thus both the average longitudinal shortening and the rotation of the end-plates 
could be determined. 


Observations through a theodolite. on a vertical line drawn on the surface of a 
cylinder, showed that the line remained vertical up to and beyond failure, thus 
indicating the absence of any torsional effects. 


7. Test Results 


The main test results have already been given in Table I. Further details are 
given in Table II. 
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TABLE II 
Size of Buckles (in.) 
Specimen | Pressure Buckling Load (circumferential x Remarks 
(Ib./in.2) (tons) longitudinal) 
1 0 7-9 | 13 x 13 Buckles at 7-9 tons; second load- 
ing buckles at 4-6 tons, third at 
| 4- 2 tons. 
2 i; O 5-9 13 x 13 Buckles could be heard forming 
at 2-8 and 5-4 tons. Load fell off 
at 5-9 tons. 
3 28-4-31-0 28-8 rDX 2 Plastic bulge formed slowly. 
Load fell off sharply at 28-8 tons. 

4 0 5-2 3x13 Buckles could be heard forming 
at 2-7, 4-6, 5-2 tons. Load fell off 
at 5-2 tons. 

5 20-7-24-6 25-4 9x 3 Plastic buckles formed slowly. 

i Load fell off sharply at 25-4 and 

6 | §8-1-12-0 15-8 9x 4 15-8 tons respectively. 

7 0 8-6 13x13 Buckles at 7-1, 8-1, 8-6 tons. 

Load fell off at 8-6 tons. 
8 0 6-5 13 x 13 Buckles at 3-6, 5-8, 6°5 tons. 


Load fell off at 6-5 tons. 


At a load of 7-9 tons, Specimen | buckled suddenly and noisily, with a marked 
fall off in the applied load. With the load steady at this lower level, existing buckles 
disappeared and new ones formed at other locations. On removal of the applied 
load the buckles apparently disappeared completely, but a close inspection revealed 
some permanent deformation at points corresponding to the corners of the diamond 
buckles. As shown in Table II, on applying the load a second and third time, lower 
values were obtained for the buckling load. 


The modes of failure for the other unpressurised tests (Specimens 2 and 4) were 
generally similar to Specimen 1. However, in both these later tests, isolated buckles 
were heard and seen to form at very low applied loads. The formation of these 
premature buckles did not seem to affect the load-carrying capacity of the specimens 
to any marked extent, so that loading was continued until a sharp fall-off in load, 
accompanied by an extensive buckle pattern, occurred. Details of the loads at which 
buckles formed are given in Table II, and the load/deflection curves for Specimens 
| and 2 are shown in Figs. 4 and 5. These figures show that the load/deflection 
characteristics of Specimens | and 2 were very different from each other. 


It has already been mentioned that the buckles in Specimen 1 substantially 
disappeared when the applied load was removed; in the case of Specimen 2 the 
buckles largely remained when the applied load was released. Some time later, how- 
ever, when the end-plates were unbolted, the buckles in Specimen 2 suddenly 
“jumped out,” leaving only small permanent deformations corresponding to the 
corners of the diamond buckles. Specimen 4 was given a considerable axial 
deflection after failure and the diamond buckle pattern was permanently formed in 
the walls of the cylinder. 


156 The Aeronautical Quarterly 


4 4 
i 


COMPRESSION TESTS ON THIN-WALLED CYLINDERS 


1 
/ THEORY c cn 
/] 
/ 
LT BO TILT 'AC 
O 
2'6 
2 ae | 
O23 Q } | 
{ 
ig ; 
2 ia ai 
, | | My, | 
a 
| | 
| 
i TILT x10 73 
10 20 3 a0 50 60 70 8 9 
DEFLECTION (in x10) 
FiGuRE 4. Specimen 1, unpressurised, 6 ft. 
7 
/ THEORY 
4 
7 FAILED 
9 
B c 
4 
y 5 j "BUCKLES 
BUCKLES | 
} 
} \ 
2 + 
| 
2,4 6 6 0.2 4 6 B 
TILTUN. x10 ~ 
10 20 30 40 50 60 70 60 


DEFLECTION (IN.x 
Ficure 5. Specimen 2, unpressurised, 6 ft. 
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In contrast to the unpressurised tests. the buckles which occurred during the 
pressurised tests (Specimens 3, 5 and 6) formed slowly. The first evidence of the 
formation of buckles came in a greater tendency for the internal pressure to increase 
due to decreasing volume; the rate at which water was being bled from the system 
had to be increased to maintain a constant internal pressure. At this stage the 
buckles could be seen to be increasing their depth in the radial direction, although 
their plan size remained sensibly constant. When the buckling loads recorded in 
Table II had been attained, a sharp falling off in applied load accompanied by a rise 
in internal pressure was observed. The buckles which formed in these pressurised 
tests remained almost unchanged after the removal of both the applied load and the 
end-plates. 

Specimens 7 and 8 were tested in an unpressurised condition by applying the 
compressive load through the parallel platen device. The general behaviour and 
mode of failure of these specimens was very similar to that already described for the 
unpressurised specimens loaded through a ball joint. The load/deflection curves 
for Specimen 7, taken up to failure, are shown in Fig. 6, and Fig. 7 shows the load/ 
deflection curves for Specimens 7 and 8, taken beyond failure up to a maximum 
longitudinal deflection of 10 in. The appearance of the specimens at large 
longitudinal deflection may be seen in the photograph of Fig. 2. 

It had been expected that the parallel platen device would prevent all rotation 
of the end faces of Specimens 7 and 8; however, this proved not to be the case. The 
end faces of these specimens did rotate, as the applied load was increased, resulting 
in the specimens being eccentrically loaded as the point of contact between the 
parallel platen and the concrete backing slab shifted from the axis of the cylinders. 
This can be seen in the photograph of Fig. 2, and is illustrated diagrammatically 
in Fig. 7. 

It has been mentioned earlier that measurements were made of the cylinder wall 
irregularities; a typical sample of the results obtained from such measurements is 
shown in Fig. 8. These measurements were taken at one-inch intervals along 
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Ficure 8. Cylinder wall irregularities (in inches). 
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selected longitudinal generators (A and B), using a 1/1000 in. dial gauge mounted on 
a slide which moved along a rigid steel straight edge. The results shown in Fig. 8 
indicate that, before test, Specimen 2 had maximum wall irregularities of + 0-020 in.; 
after buckling under zero internal pressure and removal of load, this figure could be 
+0-°130 in. The results shown for Specimen 3 indicate a very significant reduction 
in wall irregularities due to the application of internal pressure. 


The results of subsidiary tests carried out to determine material properties are 
recorded in Appendix I. 


8. Theory 


The small deflection theory put forward by Southwell”, gives a value for f...’ of 
0-605 (for a material having a Poisson’s ratio value of 0-3). Leggett, and others, 
confirmed the value for f,,,’ of 0-605 for small deflections, but showed that the stress 
after buckling, when the deflections are large, corresponds to a value of f.,,’ much 
lower than 0:605. Lo, extending the large deflection theory to include the effects of 
internal pressure, obtained a theoretical relationship connecting f,,’ and p’. This 
theoretical relationship is plotted in Fig. 9, together with the experimental results; 
the two p’ values shown for each test correspond to the pressure at the top and 
bottom of the cylinder. On this analysis f.,,’=0-376 and f.,.’ =0-605 for p> 0-17. It 
will be seen that a line drawn through the experimental points would be approxi- 
mately parallel to the theoretical line but offset some distance from it. This kind of 
observation led Fung and Sechler“ to suggest that Lo’s result should be plotted in 
the form (f..’—f..:) against p’, and this curve is shown in Fig. 10. 
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The experimental points on Fig. 10 were plotted using the relationship 
(f Fore = Ber es 0-22 (1) 


The value of f.,,,’=0:22 was derived from the following empirical relationship 
obtained by Kanemitsu and Nojima“: 


fers’ =(R/t) [9 (t/ RY + 0°16 i «. 

This gives a value f.,,’=0°22 for both the 6 ft. and 9 ft. cylinders. By using 

relation (1), it is apparent that better agreement between experimental and predicted 
results is shown in Fig. 10 than in Fig. 9. 


9. Discussion of Results and Conclusions 


The results of the five unpressurised tests show a considerable variation in the 
buckling stress, the experimental values of f.,,’ being between 0:14 and 0°24, with an 
average value of 0:19; these compare with a value of 0:22 from Kanemitsu and 
Nojima’s empirical relationship (equation (2)). 


There is some evidence suggesting that the scatter of the test results was due to 
the magnitude of the cylinder wall irregularities. This experimental evidence is 
summarised as follows : 


(a) Table III shows the average wall irregularity obtained from measurements 
of at least 100 points on the surface of each cylinder. It will be seen that high wall 
irregularity is associated with relatively low failing loads. The actual readings from 
Which the average wall irregularities were derived are shown in Appendix II. 
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TABLE Ill 
| Average Wall Irregularity Failing Load 
Specimen End Condition (1/1000 in.) (tons) 
ball 
4 ball 33 
7 platen 7 rds 
8 platen 13 6°5 


(b) Test results from Specimen | show that for successive loadings the failing 
load was 7:9, 4:6, 4:2 tons. This indicates that each successive loading increased 
the wall irregularities and hence decreased the subsequent failing load. 


(c) Figures 4 and 5 illustrate the large variation in load/deflection 
characteristics between Specimens | and 2. A relatively small deflection per unit 
load is seen to be associated with a relatively high failing load. Large deflections 
will occur where wall irregularities are high, as axial shortening due to bowing of 
the cylinder walls will be superimposed onto the shortening due to direct axial strain. 


The parallel platen was imperfect in so far as it allowed some relative rotation 
of the end faces of the specimens to take place, resulting in an eccentricity of loading 
tending to restore the end faces to their parallel positions. The following results, 
however, show an increase in strength of the parallel platen tests over the ball-ended 
tests : — 


average failing load, three ball-ended tests 6:3 tons, 


average failing load, two parallel-platen tests = 7:6 tons. 


By pressurising the cylinders, the wall irregularities were reduced and, in 
addition, the buckle pattern at failure was considerably changed. The results show 
that, as well as the pre-tensioning effect of pressurisation, there was also a very 
considerable increase in buckling stress caused by the applied internal pressure. 
Fig. 10 shows that, for internal pressures corresponding to p’ > 0:5, large deflection 
theory accurately predicts the failing load of the specimens; the agreement for 
Specimen 6 where p’ < 0:5 is not very satisfactory. 


The size of the buckles that arose in the unpressurised tests is of interest. These 
had diagonal dimensions (i.e. semi-wave lengths) of about 13 in., whereas classical 
linear theory indicates a value of about 2:7 in. Thus the experimental buckles were 
nearly five times this theoretical size. This result is qualitatively similar to that 
obtained by Robertson, but in his tests the ratio was only about two. The change 
in buckle size and shape due to pressurisation is generally similar to that observed 
by Fung and Sechler™. 
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Appendix I 


The following material properties and dimensions were measured from Specimen | 


Diameter = 36 in. 

Wall thickness =0-035 in. 

Ultimate tensile stress = 30-9 tons/in.’ 

0-1 per cent proof stress =20-4 tons/in.? 

Modulus of elasticity = 4,670 tons/in.2=10-5 x 10° Ib./in.? 
Cross-sectional area = 3-96 in.” 


Appendix II 


The actual readings from which the average wall irregularities quoted in the paper 
were derived are given in Table IV. Their derivation is illustrated by Fig. 11. 


3 INS SIN. 
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| \ Figure 11. Cylinder wall 
t | irregularities. 
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TABLE IV 
AVERAGE WALL IRREGULARITIES 
+ Outward — Inward 
Specimen 2 Specimen 4 Specimen 7 Specimen 8 
(6 ft. long) (9 ft. long) (6 ft. long) | (8 ft. 6 in. long) 
Line Row Row Row | Row 
A B Cc A B A B B 
5 0 0 0 0 0 0 0 0 0 | 0 0 0 
10 -10 ~3 +9 -1 +1 +11] 
12 — 8 3 13 —~2 +4 5 
13 +12 44 44 
14 —-10 — 4 13 -44i4 8 
15 —-9 — 3 13 -5 +5 8 
16 -9 —1 +14 -8 +7 8 
17 —-8 +1 +13 -8 +7 8 | 
18 +4 +14 +8 +9) 
19 -7 +4 +415 -8 +9 +10) 
20 +8 +16; — 1 6 1; -9 10 10; —23 -—21 
21 —-6 +7 +14 9 4|-9 10 10; -24 -21 -—9 
22 — 4 +14 17 5 9 5|— 8 10 10; —25 -19 — 6 
23 3 +13 +18; +13 146 +10; - 9 10 13 | —24 -17 -—4 
24 5 +16 +14; +15 +19 +8/-—9 4410 14; —26 —12 2 
25 |} 5 +12 16 18 22 12; -9 +10 14; -—28 -10 — 2 
26 — 6 16 18 | +18 26 3;-9 +9 12} -26 —7 0 
27 — 6 +17 +23 +30 +17/-8 +6 -2 +1 
28 6 +11 +12 +27 +31 +20) ~—10 6 10; -26 —3 +1 
29 7 IS +28 +33 +23) -10 + 6 10; +1 +3 
30 6 +12 +18 29 +35 +22|;-9 +6 +8] —% 0 +4 
31 — 7 +13 18} +31 +37 +24| -10 +7 +6] -25 +6 +8 
32 - 9 15 13 | +38 +35 +26; -11 +4 +6+9 
33 +11 +41 +35 +26] -11 +4 + 5| -22 +8 +11 
34 | +14 +13 +45 +435 31; -11 +2 +3] —25 8 11 
35 14 12; +50 +34 +3 +2] -18 11 13 
36 —19 +13 14' +35 +43 435/-9 + 3 2] -18 13 14 
37 +30 +43 +39 | -16 +14 +15 
38 24 +49 +38 —16 13 13 
39 17 52 +35 —15 13 13 
40 10 57 +34 | —13 1] +16 
41 +11 58 +35 |} —11 Il +15 
42 20 55 +34 —12 +10 +14 
43 24 55 +34 ~10 +49 15 
44 30 +54 +34 -— 7 10 11 
45 23 52 33 — 6 8 8 
46 27 49 33 
47 30 +45 30 — 4 3 +2 
48 31 +46 31 
49 +-34 44 32 
50 +28 +43 +32 
51 26 +38 29 
$2 23 +35 29 
53 24 +40 25 
54 21 +28 21 
55 17 +23 17 
56 + 9 +20 +16 
57 + 4 +14 8 
58 + 4 13 5 
59 6 +9 4 
68 0 0 0 0 0 0 0 0 0 0 0 0 
Average 11/1000 Average 33/1000 Average =7/1000 | Average = 13/1000 
in. in. in. in. 
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SUMMARY: Previous theoretical work on mass transfer cooling is reviewed and it 
is shown that this may be complemented by similar solutions that occur when 
the velocity outside a two-dimensional boundary layer varies as some power of 
the distance from the front stagnation point. The case of stagnation point flow 
with constant wall temperature is investigated in some detail, under the assump- 
tion that the temperature differences are everywhere small compared with the 
absolute temperature. Calculations on an analogue computer, supplemented by 
an investigation of the asymptotic behaviour, are used to determine the 
boundary layer development and heat transfer rates when the coolant is 
hydrogen, helium, steam or carbon dioxide. It is found that, on a mass flow 
basis, hydrogen reduces the heat transfer rate most and that steam is the next 
most effective of the substances investigated. 


l. Introduction 


The interest in mass transfer cooling has increased recently because it is a 
promising method of reducing the surface temperatures of high speed aircraft, 
missiles and satellites re-entering the earth’s atmosphere. 


In mass transfer cooling a foreign gas is emitted from the surface of the body 
and diffuses through the boundary layer, thus affecting its development and the heat 
transfer rate. The emission can arise from the evaporation of a liquid, the 
sublimation of the surface material or from the injection of a foreign gas through 
holes in the surface. 


Extensive theoretical work on this topic has been reported for the case when the 
molecular weights of the emitted gas and that in the external stream are the same, 
and this is summarised in Ref. 1. Analytical investigations of the more general 
problem appear to have been confined to the case when the Mach number in the 
external stream is constant* and the emitted gas is a light one, either hydrogen or 
helium (e.g. see Ref. 3). 


In the present paper it is shown that this work may be complemented by 
similar solutions that occur when the velocity outside a two-dimensional boundary 
layer varies as some power of the distance from the front stagnation point. The case 


*When this work was nearly completed the author's attention was drawn to Ref. 2, which 
deals with the effect of helium injection on the heat transfer near a stagnation point in 
axisymmetric flow. 
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of stagnation point flow with constant wall temperature is then investigated in some 
detail, under the assumption that the. temperature differences are everywhere small 
compared with the absolute temperature. 


NOTATION 
A constant (see equation (49)) 
(see equation (57)) 
B constant (see equation (49)) 
b=mMS,, (see equation (58)) 
C constant (see equations (17) and (49)) 
C,,C.,C,,C, constants (see equations (44) and (48)) 
C, specific heat at constant pressure 
constant (see equation (51)) 
D,, coefficient of binary diffusion 
E constant (see equation (54)) 
F constant (see equation (57)) 


Poo¥ 


f non-dimensional stream function, rene 
xX?’ 


1 
Seo an 
G=g-n+ 8 
h=2+mM (see Section 6.5) 
H, stagnation enthalpy in external stream 
k coefficient of thermal conductivity 
1=0, —%p w (see Section 6.5) 
m=m,/m, 
m,,m, molecular weight of foreign gas and air respectively 


M’ mass flow per unit area of surface per unit time 


1/2 
M_ mass flow parameter, = M’ = | 
n constant (see equation (17)) 
(Nu) Nusselt number, Ox 
k.(T.,-Tw) 
pressure 


(Pr) Prandtl number 


Q heat transfer per unit area of surface per unit time 
r see equation (62) 
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R_ (in Sections 2 and 3 only) universal gas constant 
| R_ (elsewhere) Reynolds number. = uy 
S Schmidt number, =v/D 
T temperature 
u component of fluid velocity in direction Ox 
v component of fluid velocity in direction Oy 
w mass fraction of foreign gas, =p, /p 
x distance measured along surface from stagnation point 
y distance measured normal to surface from wall 
& constant (see equation (39)) 
» non-dimensional distance normal to surface, =) | 
UX 
6=9,-1 


v 


MASS 


coefficient of viscosity 


kinematic viscosity 


- = —8 (see equation (39)) 


fluid density 


D, 


TRANSFER 


OOLING 


Cz 


2a0° 


respectively 


ky, k 


k,, 


respectively 


Or. Pus %, | Lec p respectively 
fe 


stream function 


Suffixes 

| denotes foreign component in binary mixture 
2 denotes primary component in binary mixture (air) 
co denotes value at outer edge of the boundary layer 


w_ denotes values at wall. 


2. Boundary Layer Equations 


The development of a two-dimensional boundary layer, on a surface through 
which a foreign gas is being ejected, is considered. When conditions are steady the 
equations which determine the fluid motion, temperature and concentration of 
foreign gas are equations (1)-(6) which follow.* 

*Equations (1) to (5) are identical with those given on pp. 456-7 of Ref. 4, except that the 
sign of the second to last term in the energy equation is misprinted in that reference 
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Continuity : 
0 0 
Diffusion: 
ow Ow ow 
pu ox f oy dy | Dir 4 (2) 
Momentum: 
Ou Ou op. Ou 
) + pl = = . . 3 
and dp/dy=0. ; (4) 
Energy: 
oT oT op 0 
dx * dy 0) dy 
ow OT 
+ 5 
+ (C,, Cy2) ey oy iz) (5) 
State : 


p Mm, m, 


In these equations x and y are distances measured around the surface from the 
front stagnation point and normal to the surface respectively, and u and v are the 
velocity components in these directions. p, ,C,, T and k are the density, coefficient 
of viscosity, specific heat at constant pressure, temperature and coefficient of thermal 
conductivity of the gas mixture respectively; w is the mass fraction of the foreign 
gas, D,, the coefficient of diffusion, p the pressure, C,, and C,, the specific heat of 
the foreign gas and air respectively, m, and m, their molecular weights and R is (in 
this section and in Section 3) the universal gas constant. 


In general the physical properties D,,, u, C, and k will depend on both w and 7, 
and C,, and C,, on T alone. Subsidiary equations defining these dependences are 
needed to render the six main equations soluble for p, u, v, w, p and T. 


It is noted that the diffusion equation as stated does not involve the assumption 
that w is small and that, in the energy equation, account is taken of the energy 
transfer that arises from diffusion processes. 


3. Boundary Conditions 


When the foreign gas is ejected through holes in the surface, the mass rate per 
unit area, M’, at which it is transported by both convection and diffusion, may be 
considered to be prescribed. 
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Thus, the boundary conditions at y=0 are as follows: 


(ii) The rate per unit area at which air is transported is zero. Thus 


Ow 


pv (1 —w)+ D,.p =0. (8) 


Addition of equations (7) and (8) gives 


which, when substituted into equation (7), gives 


~ 


=M’(1-w). 


(iii) The tangential velocity component is zero, so that 


(iv) The wall temperature is prescribed, so that 


where suffix “w’’ denotes values at the wall y=0. 


It will be assumed that the heat transfer rate per unit area of surface, Q, is 
given by 


When the mass transfer arises from evaporation or sublimation, M’ is not known 
and an extra boundary condition relating it to the conditions at the wall must be 
satisfied. Thus the solutions for these cases may be obtained from the ones when 
M’ is prescribed, which are considered later, by selecting those solutions which 
satisfy this additional condition. 

For all cases the boundary conditions at the outer edge of the boundary layer 
are the same and are, as y —> ©, 


u— Us, 
and T—>T, 
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Suffix 0© denotes values at the outer edge of the boundary layer. u, and T, 
must satisfy Bernoulli’s equation, which may be written 


where H, is the stagnation enthalpy in the external stream. 
At the outer edge of the boundary layer, w = 0, so that equation (6) gives 


Px _ ‘ 
Po m, 


Combining this with equation (6) gives, in virtue of equation (4), 


oo, [1+ (m—-1)w]=1, — | 
where ¢,=p/p.~ and m=m,/m,. 


4. Similar Solutions 


It would appear that similar solutions of the compressible boundary layer 
equations may only be obtained under one of the following three conditions : 


(i) That the Mach number of the external stream is constant. 


(ii) That the property values of the gas are assumed to obey certain 
simple laws: e.g. «OC T and the Prandtl number, (Pr), is unity 
(see Ref. 5). 


(iii) That the Mach number of the external stream is small compared 
to unity. 

If condition (i) is satisfied, then similar solutions may also be found to 
equations (1) to (6) of Section 2, which include the effects of mass transfer“. 

In the present paper it will be assumed that condition (iii) is satisfied, that the ' 
wall temperature is constant and that the velocity outside the boundary layer is 
given by 


where C and v7 are constants. 


Equation (15) may now be replaced by 


Poo + = constant (18) 
and 7.,=constant. 


Equation (1) is satisfied by introducing a stream function, v, such that 


Px oy / dy ) 


and pv= --p,, { ° (IY) 
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Let (20) 
and f- Po (21) 


and suppose that f is a function of » only, so that 


_ Ux af 
dn | 
(22) 
df n+l 
and dh ‘|. | 


where the Reynolds number R = pu. X/ 


The assumption that the Mach number is small enables the first and last terms 
on the right-hand side of equation (5) to be neglected so that, using equations (17) to 
(22), equations (2), (3) and (5) become, respectively, 


n+1, dw dw 
n (df\* (n+1) ,d {1 d d {1 
24 
2 dn dn dn\¢, an (24) 


0 iw d 
and (Pr). (n+ 1) d (Pr), dw doy 


? fo, dn dn Dx dn = 9 12 %, (¢ 12 1) dn dn ( ) 


where Discs Hoos 


(Pr). pos Rap $c =C, / and = Cu 


The boundary conditions, equations (9) to (12) and (14), become, at »=0, 


M (1-w)= S. (26) 
df /dn=0 
and 
where M = M’ = | ig is the mass flow parameter, and, as » —> ©0, 
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1 
and or— 


For the four equations (16) and (23) to (25) to be soluble for f, w, 9, and 9r. 
subsidiary equations are needed which define the dependence of the property values 
9p and on w and or. 


5. Stagnation Point Flows for Small Temperature Differences 


For the present it will be assumed that the temperature differences are every- 
where small compared with the absolute temperature, so that equation (16) 
simplifies to 


(28) 
It will also be assumed that 


which is valid near a stagnation point where the heat transfer problem is particularly 
severe. 
By equations (28) and (29), equation (23) becomes 


f do, 


it being assumed that #,,, is unity, an assumption that is well satisfied for all gases 
since D,, is nearly independent of concentration. 


It is convenient to define a new function ¢ by the relation 


1 d9, 
(31) 
ay 
so that equation (30) becomes 
do, 1 d*9, 
32 
dn 
Differentiation of equation (31) gives, using equation (30), 
dg 1 df 
; 33 
dn dy 


From equations (29), (31) and (33), it follows that equations (22), (24) and (25) 
become equations (34)-(36) which follow. 
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ag 
= Ugg a 
(34) 
de\? d*g 8) 
in) dy? \ dy %p (* dn? 65) 
d*¢, doy | doy do,\ (Pr), (1—%e12) 
The boundary conditions (26) and (27) become, at  =0, 
>= —mM 
dg|dn=0 | (37) 
dr — | | 
1 (38) 
| 


The assumption that the temperature differences are small compared to the 
absolute temperature enables the dependence of the property values ¢,, 9, and 9%, 
on ¢7 to be ignored. Thus equations (32) and (35) do not involve ¢; and may be 
solved for g and ¢,. Once their solution is known, equation (36) may be solved 
for 


6. Solution of Diffusion and Momentum Equations 


The solutions of equations (32) and (35) were investigated by using an 
approximate method and also by calculations on an analogue computer. Since each 
approach involves a knowledge of the asymptotic behaviour this will be 
investigated first. 


6.1. ASYMPTOTIC BEHAVIOUR 


G, € and 4 are defined by the relations 


f=r-8 | (39) 
and 9¢,=1+4 


where 4 is the constant which is such that, for 1 large, G is small (see the first of 
boundary conditions (38)). 
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Substitution from equation (39) into equations (32) and (35) gives, respectively, 


dé 


dé = —§,,(€+G) . (40) 
and 
dG dG\’ | dé d’G 
(FF) Cae ~ dé dé? 
146° 146 dé d& 
In equation (40), G is neglected compared to €, so that the equation becomes 
dé A? 


In equation (41) it is supposed that G, 4, (1 —¢,,) and their derivatives are small and 
only terms of the first order are retained. Thus the equation becomes 


dé -$_£7/2 
where C, is a constant. 
2 
2 
where erfc x= | 
x 


and C, is a constant, which is not taken to be zero at this stage, although this implies 
that the second of boundary conditions (38) is not satisfied. 


From equation (45), equation (43) can be written 


@G .@G Ste 
de + qe - C: etfc +C;. (46) 
Differentiating this equation twice gives 
dé +€ -C,S,£ e . (47) 
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It follows that, when s,, = |, 
C S_ -£2/2 


+ constant x e 
dé -S,. 


which may be integrated three tines to give 


dG eu 
dé Uy 
5.7 (-S,)' erfc 3173 +C,i +C,(€ +1), . (48) 


where /"erfc.x | i”  erfc dt, 


x 


with i erfex=erfc x, and C, and C, are constants. If the first of boundary 
conditions (38) is satisfied, then C, and C, are zero. 


When S,=1, 


d‘G 


dé' 2 


+constant x e 


so that the asymptotic behaviour is not very different from that when S, = | and 
will not be considered further herein. 


6.2. PROPERTY VALUES 


Equations (32) and (35) were solved on the A.R.L. analogue computer for a 
range of values of M when the ejected gas was hydrogen, helium, steam or carbon 
dioxide; but, before considering these solutions, the property values used in the 
calculations, and also in the subsequent solution of equation (36), will be discussed. 
They are given in Table I. 


The values of »,, k, and D,, for the gases under consideration were obtained 
from Refs. 6 and 7 for a pressure of one atmosphere and a temperature of 0°C, 
except for steam*. These values were used to give values of %,,., %,, and S, for 
each gas, where 

Puyo = My] Me and 12 K, | Re: 


*For steam the property Vafues, and also the reference property values of air, were evaluated 
at 100°C, except that, in evaluating Purar My WAS erroneously taken as the coefficient of viscosity 


of air at O°C. This gave a value of 0°743 for u,/u,. According to Ref. 8 (Table 2V-3) the 
value of this ratio with both quantities evaluated at 100°C is 0-605. It is believed that the 


effects of this error are small. 
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TABLE I 
PROPERTY VALUES 

Property | Hydrogen Helium Steam Carbon Dioxide 

m 14-48 7-328 1-608 0-658 
Ss 0-216 0-254 0-600 0-957 
| 0-4927 1-077 0-743 0-802 
A 0-4145 1:0 0-3203 1-3817 
B | 1-1708 0 0-6797 ~0-3817 
Cc ~0-5845 0 0 0 

7-13 6-23 0-7663 0-614 
Gers 14-13 5-17 2-005 0-829 
(Pr), 0-72 0-72 0-72 0-72 

Range of 0 to 0-4 0 to 0-6 0 to 1-0 0 to 1-6 


Reference 6 also gives the variation with concentration of the coefficient of 
viscosity of a mixture of hydrogen and air and of a mixture of oxygen and carbon 
dioxide. It was found that, in the former case, the values could be closely 
approximated to by 


the relationship between ¢, and w being given by equation (28). The values of A, B 
and C are given in Table I. In the latter case a linear approximation was found to 
be satisfactory. In view of this, and because of the dearth of other experimental 
results, it was assumed that a linear formula could also be used for mixtures of air 
and carbon dioxide and of air and steam. The values of the constants A and B were 
chosen to give the correct values of ¢, for pure air and pure foreign gas and are 
included in Table I. The coefficients of viscosity of helium and air were considered 
to be sufficiently close for it to be assumed that the viscosity of a mixture of them is 
independent of concentration. 


For mixtures of hydrogen and air and of helium and air, the formulae given in 
Refs. 3 and 2 were used respectively to give the dependence of ¢, on ¢,. For 
mixtures of steam and air and of carbon dioxide and air, the variation of ¢;, with @, is 
much smaller and the simple formula 


(x,2.— 
Dy l 
m-1 
was used. This formula gives the correct values of ¢, for pure air and pure foreign 
gas, when ¢,=1/m (see equation (28)). 


6.3. CALCULATIONS ON THE ANALOGUE COMPUTER 

The calculations on the analogue computer were commenced at »=0 with 
prescribed starting values for g, dg/dn, ¢,.d*2/dy?, ¢, and do,/dn. The starting values 
of the first two functions were obtained from the first and third of equations (37). 
The other starting values were determined as follows. 
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(i) Starting values of , and ¢, d?¢/dn? were guessed and the second 
of equations (37) was solved for d¢,/ dv. 


(ii) The solution curve of dg/dy for these starting values was then 
obtained. In general it was found that, for large values of », 
dg/dn tended to either +00 at a rate which appeared to be 
consistent with the last term in equation (48). The starting 
value of ¢, d*g/dn? was then adjusted until the divergence of 
dg/dn towards +°0 was delayed to as large a value of » as 
possible. 

(iii) The solution curve of ¢, was then obtained. In general it was 
found that, for » large, ¢, tended to a constant other than one in 
accordance with equation (45). 


(iv) The starting value of ¢, was then adjusted to a value that was 
estimated to be such that, for » large, #, would tend to one. A 
new starting value of do,/dy was obtained from the second of 
equations (37). 


Procedures (ii), (iii) and (iv) were then repeated in turn until a solution was 
obtained that satisfied the two boundary conditions at »=°°O with acceptable 
accuracy. 


Typical solution curves that were obtained in the final stages of calculation for 
the case of helium injection with M=0°3 are shown in Fig. 1. The figure shows the 
effect of making a small change to the starting value of ©, d*g/dn?, keeping the 
starting values of the other functions the same. 


25 


0-75 


0°50 


025 


ce) 


Figure 1. Typical solution curves obtained on an analogue computer showing the effect of a 


small change in the starting value of ¢, d?g/dy?, keeping the starting values of ¢,, d?,/dn, g 
and dg/dy the same. Helium injection, M=0°3. 


177 


| | | 

| | | | | 

| | | | | 


D. BURLEY 


The values of 9, (and d¢,/d») are not affected and it is seen that the boundary 
condition at 7 =°° on ¢, is satisfied accurately (to within approximately 0-001). 


The values of dg/dv are affected little for 1, < 3, but for 1 > 3 the two curves 
diverge, one tending to +00 and the other to — 00. Thus the correct starting value 
of », d*g/dn? can be determined very accurately, but the correct values for dg/dn 
for » > 3 are left in doubt. It was found that these values could be determined with 
acceptable accuracy by using the knowledge of the asymptotic behaviour of the 
solution (see Section 6.1). The values of dg/d1 obtained in this way are shown by 
the dotted curve in Fig. 1. 


This instability in the values of dg/dn for 4 large is typical. It became more 
pronounced as M was increased and imposed an upper limit on the value of M for 
which solutions could be obtained for a particular gas. 


The results obtained on the analogue computer, corrected for » large, as already 
explained, are given in Figs. 2, 3 and 4. 
6.4. Exact SOLUTIONS FOR SPECIAL CASES 
Equations (32) and (35) may be solved exactly for two special cases : — 
(i) for M small 
and (ii) for m=1 and M arbitrary. 
In both cases @, and @, may be taken as unity in equation (35) so that it becomes 


5 
dn) 8 dn? i+ dy’ (50) 


which is the familiar equation for the non-dimensional stream function for 
stagnation-point boundary-layer flow, solutions of which are available. 


For M small, equation (32) is solved for ¢, to give 


" 
[ exp | e@ar) dé, 


where D is a constant. 


From this equation it follows that 


[exp [ at) ac 
| exp | (t) at) dé 


0 0 
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FiGurE 2. Boundary layer velocity profiles for injection of various gases. 
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FIGURE 3. Boundary layer density profiles for injection of various gases. 
~ — Analogue solutions. 
—— --—— Exact solution for limit M>0. 
—-—--—— Approximate solutions (see Section 6.5). 
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| | m =14: 
| J 0-216 
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0-2 O06 M 08 2 4 
Ficure 4. Variation of concentration at wall with mass injection parameter 
for injection of various gases. 
O-—— Analogue solutions. 
Exact solutions for m= 1. 


a Approximate theory for helium injection 
(see Section 6.5). 


The values of g for zero fluid injection were obtained from Ref. 9, and equation (52) 
was evaluated for the values of S,, appropriate to the gases being considered. The 
tesults are included in Fig. 3. 
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For m=1 and M arbitrary, equation (28) is used to transform equation (32) into 
one for w, the limit »—> | being taken after the transformation. Thus 


dw 
dy? S08 dn ( 
the solution of which is 
w=E | exp | (t) dt) dé, 


where E is a constant. Boundary conditions (37) yield, at »=0, 


g=-M 
dw /dn= —MS,, (1 (S5) 
and dg/dy=0 | 
From the second of these, and from equation (54), it follows that 
MS,, | exp ( | dt) dé 
We = = . (56) 


~ é 


| exp ( g(t) dt) dé 


0 0 


1+MS, 


where w,, is the concentration at the wall. 


The solutions of equation (50), subject to boundary conditions (55), are given in 
Ref. 10. These, together with the same values of S,, as used in the analogue calcula- 
tions, were substituted into equation (56) to give the dependence of w, on M. The 
results are included in Fig. 4. 


6.5. APPROXIMATE METHOD OF SOLUTION WHEN ¢@, = 1 


An approximate method of solving equations (32) and (35) when 4, = 1 was also 
investigated. It consists of matching, at a point within the boundary layer, the 
asymptotic solution derived in Section 6.1 with a solution which is valid near the 
wall. For » small, put 


g=—mM+h 
and ¢,=¢pwtl. 


It is supposed that /, / and their derivatives are small and only terms which are 
of the first order in these quantities are retained in equations (32) and (35). 
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Thus these equations become, respectively, 


d?l dl 
dr? mMS,. dn 
and mM bye 1=0. 
dn dr? 


The solution of these equations which satisfies boundary conditions (37) is 
e=—mM+h mM + F (1 + an —e™)+ (57) 


where a= Mo,w, b=mMS,, and F is a constant 


From equations (45) and (48), the asymptotic solution which satisfies boundary 
conditions (38) is 


P i Sxo : 5 


and g=¢- .C,Perfc -2""C, Perle =, . (60) 

It is supposed that this solution is valid for €>0, ie. for 7» >4 (see 

equation (39)), and that the solution given by equations (57) and (58) is valid for 

» <6. In matching the solutions at »=4, there are five disposable constants, 6, 

C,,C, and F. These were determined by requiring ©,, d¢,,/dn, g, dg/dn and d*g/dn? 
to be continuous at = 64. 

The boundary layer velocity and density profiles were calculated for helium 

injection using the property values given in Table I. The results are compared with 

the corresponding analogue ones in Figs. 2, 3 and 4. The comparison provides a 


check on the accuracy of the approximate theory, as in both sets of calculations 
was taken to be unity. 


%7. Solution of Energy Equation 


The energy equation (36) can be written in the form 


dor 
-(), ; (61) 
1 | do, 1 do, (Pr), (1 -%ei2) 4%, (62) 
where » (n)= + (Pr),. (2 + wy dn 
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Ficure 5. Boundary layer temperature profiles for injection of various gases. 


184 The Aeronautical Quarterly 


| //\/ | | 
TTT 
0-4 ‘ow ff | 
YZ | | | | | | 
1 
| 
0-054 
0-6 xa 
I 
1 
+0 | 
Yj 
0-4 7 
| 
0-2}- Ad | 
| 


MASS TRANSFER COOLING 


! | 

_ 
\ 

NT AIR 

| 
0-1] H, 

H.0 

02 03 0-4 05 06 07 

M 


FiGURE 6. Variation of heat transfer rate with mass injection parameter for various gases 


The solution of equation (61) which satisfies the boundary conditions (37) and 
(38) can be expressed in the form 


exp | [ r (1) dt) dé 


(63) 
-T. 
exp r(t) dt) dé 
lhe Nusselt number, (Nu), is defined by 
Qx (64) 

and it follows from equations (13), (20) and (63) that 

(Nu) (dT dn) (65) 


T~—Tx 3 
| exp | | r (1) dt) dé 


Equations (62), (63) and (65) were evaluated numerically, using the results of 
the analogue calculations, and the resulting values are given in Figs. 5 and 6. 
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8. Discussion of Results 


To test the accuracy of the analogue computer, the values of u/u,, that were 
obtained for M=0 are compared in Fig. 7 with the exact values obtained from 
Ref. 9. The agreement is seen to be satisfactory. 


The boundary layer velocity and density profiles are shown in Figs. 2 and 3. A 
number of their features can be accounted for in terms of the asymptotic solution 
derived in Section 6.1. For each of the gases under consideration, S,, <1 (see 
Table I), so that it follows from equations (45) and (48) that, for » large, ~ 


Ux 21/2 
(66) 
(n-8) 


Thus the density and velocity profiles will be thicker the smaller the value of S,. 
This feature is apparent in the figures. 


When m > | the densities within the boundary layer will be less than the free- 
stream density. Thus, from equations (66), C, will be positive and the maximum 
value of u/u, will exceed unity. Fig. 2 shows this feature and that it is more 
pronounced the larger the values of m and M. 


Figure 4 gives the dependence of w., on M for the various gases and includes 
results derived by the method of Section 6.4 for m=1 and the same values of S, as 
used in the analogue calculations. It is seen that for the light gases it is essential to 
take the correct value of m to obtain reliable values of w,. If it is assumed that 
m= 1, then the predicted values of w, are too small, particularly for the larger values 
of M. For the heavier gases, the assumption that m=1 leads to smaller, although 
significant, errors in wy. 


Figure 4 also shows that, for all values of M, w, increases with S,, for m fixed 
(=1). This is to be expected as an increase in S.=v,/D,., corresponds to a 
decrease in the diffusion coefficient D,.... 


The results obtained for helium by the approximate theory described in Section 
6.5 are included in Figs. 2, 3 and 4. The points in the boundary layer where the 
outer and inner solutions were matched are indicated. Fig. 2 shows that the 
approximate theory gives unreliable values of u/u,, for the larger values of M. A 
possible reason for this is that the assumption in the approximate theory that u/u. 
is small over the inner portion of the boundary layer does not hold. 


The values of (p — pw)/(px,— pw) given by the approximate method are seen to be 
in closer agreement with the analogue results (Fig. 3), as are the values of wy (Fig. 4). 


The temperature profiles are shown in Fig. 5 and a noticeable feature is the 
progressive decrease in the temperature gradient at the wall as M is increased. 
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FiGuRE 7, Comparison of values of u u,., given by analogue computer with exact values. M=0. 


The values of (Nu)/R''* are given in Fig. 6, which includes results for air 
injection obtained from Ref. 10, corrected to a Prandtl number of 0°72. Two factors 
determine the value of (Nu)/R'!?. These are the thermal conductivity of the gas in 
contact with the wall and the temperature gradient there (see equation (65)). Both 
steam and carbon dioxide have lower thermal conductivities than air (see Table I) so 
that in these cases both the thermal conductivity and the temperature gradient 
decrease as M is increased. Thus, so does (Nu)/R''?. Hydrogen and helium, how- 
ever, have thermal conductivities considerably larger than air, so that for these gases 
the thermal conductivity at the wall increases steadily as M is increased. Fig. 6 
Shows that this effect outweighs the effect of decreasing temperature gradient for 
small values of M. Thus (Nu)/R'!? increases with M at first and then decreases. 


The figure shows that the ejection of hydrogen reduces the heat transfer rate the 
most, and that steam is the next most effective of the substances investigated. 


Vay 1961 187 


09 ! | 
an 
| / 
0:7; 
= | | 
| 
| 
u/u,, | 
0:5} 
| | | 


9. Possible Extension 


The effects of ionisation and dissociation of the ejected gas have not been 
considered in the present investigation. Since a very large quantity of heat is 
absorbed in these processes, it seems likely that account should be taken of them in 
determining the best substance to use for mass transfer cooling. It may prove 
practicable to extend the present method of analysis to include these effects. 
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The Elastic Stress Concentration Factors in 
Shouldered Shafts 


ParRT I: Shafts Subjected to Torsion 


I. M. ALLISON 


(Department of Civil and Municipal Engineering, University College, London) 


SUMMARY: The photoelastic “stress freezing” technique has been employed to 
evaluate the elastic stress concentration factors associated with the fillet blend 
radius in a number of shouldered shafts. The full range of practical sizes of 
blend radius and depth of shoulder has been examined. Comprehensive 
results for the shaft subjected to torsion, pure bending and axial load are given 
in this (Part I) and in two subsequent papers (Parts II and III). 


The accuracy of the graphs of stress concentration factors is better than 
six per cent. Comparison has been made with the existing theoretical and 
experimental results for each mode of loading. The results of an investigation 
into the limiting value of the stress concentration factor (for a particular 
shoulder radius) as the depth of the shoulder is increased to infinity are 

included in Part III. 


1. Introduction 


A designer is often required to choose a suitable blend radius at a change in 
the section of a cylindrical shaft which may be subjected to large static or 
alternating loads. Although other considerations, such as the notch sensitivity of 
the shaft material in the case of alternating load, may influence the design, a 
knowledge of the elastic stress concentration factors (s.c.f.) associated with different 
root radii is of considerable aid in the selection of the most suitable radius for a 
particular purpose. This applies especially in examples, such as a locating 
shoulder, where it is an advantage to use the smallest possible root radius. 


Theoretical and experimental stress concentration factors associated with a 
number of particular blend radii are available, but the information from different 
sources is conflicting and it does not cover the full range of practical dimension 
parameters. 


The investigation described in this series of three papers forms part of a 
programme initiated by the Structures Committee of the Royal Aeronautical Society 
in order to provide complete information about the peak elastic stresses associated 
with a number of different types of geometric stress raiser. In particular, the present 
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Ficure |. Shouldered shaft notation. 


work is intended to provide systematic data for the full range of shaft diameters and 
shoulder root radii which would be likely to occur in practice. Tests covering this 
range of size parameters were carried out for three different forms of applied load, 
namely, torsion, bending and axial load, and each of these will be described 
separately. Since much of the experimental procedure is common to all the forms 
of loading, a general description of the method is presented first. The problem 
involves a three-dimensional stress system and has been approached by the 
photoelastic “stress-freezing” technique. 


NorTaTION (See Fig. 1) 

D large diameter of shaft 
d_ small diameter of shaft 
r shoulder radius 
h_ shoulder depth (thus D=d+2h) 

r,o,z cylindrical polar co-ordinates 
z axis of shaft 

P,Q,R_ principal stresses 

p.q.r directions of principal stresses 
n direction of the normal to the shaft surface 
ss normal stress tangential to the surface in a direction s 


Ss cross-radial shear stress acting on the plane of the cross section 
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STRESS CONCENTRATION IN SHOULDERED SHAFTS 


The stress concentration factors are based on the nominal stress in the smaller 
diameter of the shouldered shaft. The factor K used throughout this paper is 
therefore defined as 


peak stress in the component 


maximum stress in a shaft of the smaller diameter’ 
under the same loading conditions 


2. Method of Tests 
2.1. MANUFACTURE AND LOADING OF MODELS 


The range of photoelastic models tested is shown in Table I. 


TABLE I 
h/D=0°05 Ol 0-2 0:25 03 0-475 
d/D=0°9 0-8 0-6 0:5 0-4 0-2 0-05 
1-0 x x x | 8 in. long 
2:0 x x models 
30 NP x 
40 NP x 
50 T 
6°67 m4 | 10 in. long 
75 x j models 
15-0 x | 
95-0 


Torsion test only. 
NP Photoelastic test not practicable. 


The models were machined from oversize cylindrical castings of Araldite B, 
hot-setting resin, and an annealing cycle was carried out both before and after 
rough machining. An accuracy of +0-002 in. was maintained in the dimensions 
of the finished models, and the initial casting stress was never permitted to exceed 
one fringe per inch (i.e. about 3 per cent of the peak stress). 


With the exception of the shaft having dimensions specified by the parameters 
h/r=95, h/D=0-475, all the models were made from single castings, the overall 
dimensions being chosen so that the end effects of the loading system should not 
disturb the fringe pattern in the vicinity of the shoulder. 


The minimum convenient size of shoulder root radius which may be examined 
photoelastically is limited to 0-05 in., giving, for the particular model just 
mentioned, values of D=10 in. and d=0°S in. for the larger and smaller diameters 
of the shaft respectively. It was known to be impossible to make such a model 
from one casting and still maintain the initial stress level below one fringe per inch. 
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Each of the models of this type was made from two castings, the manufacture and 
subsequent loading cycle being as described in Section 6. 


In the original programme, models having dimension parameters h/r=3, 
h/D=0-05 and h/r=4, h/D=0-05 were specified. However, in these cases the 
choice of D=1°5 in. gives values of r of 0:025 in. and 0-0188 in. respectively. It 
was decided that these shoulder radii were too small to be produced sufficiently 
accurately or to permit accurate measurements of the s.c.f., and the results for these 
models were obtained by extrapolation of the graphs. 


A special loading frame which had been used in a previous part of the 
programme” was available for loading the models, which could be subjected to 
tension, compression, bending or torsion by using different arrangements of links 
and dead weights. 


2.2. PEAK STRESS MEASUREMENT 


After freezing, slices were cut from the models and placed in an immersion 
fluid, having the same refractive index as Araldite, for examination in a standard 
photoelastic bench. The peak stress was found to lie on the surface of the shaft 
and its magnitude was obtained for all types of loading by using the method outlined 
in the following paragraph. 


The stress difference was measured at a number of points in the slice below the 
shaft surface and the peak stress was obtained by graphical extrapolation to the 
point of maximum stress on the boundary. It is necessary to adopt this extra- 
polation method because of the difficulty in locating the photoelastic “fringe” on 
the extreme boundary. The presence of a “time-edge stress,” owing to the absorp- 
tion of water in the surface of the photoelastic model, presents further difficulties 
in the accurate measurement of the surface stress. Optical measurements made 
using the Sénarmont method gave the extrapolated value of the boundary stress to 
an accuracy of 0:5 fringes per inch, while the position of points on the model was 
measured with a travelling microscope to an accuracy of +0-001 in. Measure- 
ments could not be made within less than 0:01 in. of the boundary, owing to the 
curvature of the shaft surface, and the extrapolation over this distance is 
estimated to give an accuracy in the peak stress measurement of +4 per cent. 


2.3. MEASUREMENT OF NOMINAL STRESS 


No attempt was made to measure the load applied to the model, the maximum 
stress in the smaller diameter of shaft being obtained in the uniform stress region 
lying between the disturbances caused by the shoulder and the loading pins. This 
method had two advantages in that it provided a check that the end conditions were 
not affecting the stress distribution around the concentration and it also dispensed 
with the necessity for finding the fringe value for each model. The maximum 
uniform stress was measured to an accuracy of +2 per cent, giving an overall 
accuracy of +6 per cent for the final s.c.f. 
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STRESS CONCENTRATION IN SHOULDERED SHAFTS 


Ficure 2. Slices cut from the frozen photoelastic models. 


The direction of viewing q is inclined at 45° to the plane of the slice for the torsion models 
and is perpendicular to the plane of the slice for the tension and the bending models. 


3. Torsion Tests 
3.1. STRESS PATTERN IN TORSION 


The maximum stress on any cross section of the shaft is the 5@ shear stress at 
the surface. This shear stress may be represented by principal stresses P and Q, 
lying in the tangent plane to the surface and inclined at 45° to the axis of the shaft 
such that P= —Q (Fig. 2). 


3.2. PEAK STRESS MEASUREMENT 


The peak surface stress may be obtained from any longitudinal slice containing 
a diameter of the shaft (Fig. 2). The slice was cut by hand and its surface was 
finished with fine emery paper. To measure the true principal stress difference it 
is necessary to view the slice obliquely at an angle of 45° to the normal direction. 
For convenience the slice was rotated in a model holder until OC, the normal to 
the boundary passing through the peak of the concentration, coincided with the 
microscope traverse. Observations were then made at +45° and at —45° to 
obtain (P—R) and (Q—R), the tensile and compressive principal differences, 
respectively. A similar series of observations was made for O’, the corresponding 
concentration on the opposite boundary of the slice, and the mean of the two 
tensile and two compressive results was plotted in fringes per inch against the 
distance from the shaft surface, the resulting curve being extrapolated to obtain the 
boundary stress. Taking the numerical mean of the four results eliminates the 
initial error in setting up the slice, any bending induced by the loading system and 
the time edge-stress effect. A typical extrapolation is shown in Fig. 3. 
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Ficure 3. Typical torsion test results. 
Relative retardation through an axial slice, along the line OC, 
for oblique viewing at angles of incidence of +45° and —45°. 
x View at +45° (tensile principal stress). 
@ View at —45° (compressive principal stress). 


3.3. UNIFORM STRESS MEASUREMENT 

Measurements of the relative retardation on the principal planes, which are 
inclined at + 45° and —45° to the plane of a longitudinal slice, were made at points 
along a diameter of the shaft. The results were plotted and extrapolated to the 
surface to obtain the maximum principal stress in the smaller diameter of the 
shouldered shaft. 


3.4. RESULTS 


Initially graphs of the s.c.f. K were plotted against the following parameters : - 
h/r for constant A/D, 
h/D for constant h/r, 
h/d for constant r/d, 
r/d for constant h/d, 


and d/D for constant r/d. 


All these curves must be smooth and continuous, and corresponding points on 
each curve must be consistent. The graphs were drawn to satisfy these conditions, 
while paying due regard to the experimental points, and the final curves of K 
against h/r and K against d/D are given in Figs. 4 and 5. The parameters specify- 
ing the dimensions of the “insert” model are considerably different from those of 
the other models tested. The results of this test were therefore not taken into 
consideration when determining the general shape of the curves given in Figs. 4 
and 5. 
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FiGure 4. Stress concentration factor K against //r for constant values of 4/D 


For all the shouldered shafts tests. the curves within the range of the 
experimental results have been drawn as full lines. Results for //d=0-05, outside 
the test range, were obtained by extrapolation and the curve in this region is shown 
as a broken line. 


For given values of r and d. K increases initially. reaches a maximum value and 
then decreases as D is increased. This behaviour, which cannot be explained by any 
simple argument, also occurs in pure bending and to a lesser extent in the case of 
axial load. 


For all three forms of loading and for all values of r/d, the results of the 
photoelastic tests inevitably predict an ultimate decrease in the s.c.f. as the value of 
d/D is decreased. This fall in the s.c.f. is still present even if the maximum 
estimated error band is put on the experimental points. Furthermore, were this 
trend not really present. it seems most unlikely that it should appear consistently 
throughout a relatively large number of photoelastic tests, conducted for a number 
of different forms of loading. 


For the sake of completeness the results of the insert test have been included 
in Fig. 5, and a possible shape for the complete curve for r/d=0-l is indicated by 
the chain-dotted line. It is interesting to note that the results of this test confirm 
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FiGuRE 5. Stress concentration factor K against d/D for constant values of rd. 


the trend already observed, but it should be emphasised that this curve and the 
single test on which it depends are subject to larger possibility of error than the 
family of curves in the region d/D > 0-4. 


A more complete discussion of the significance of the fall in the s.c.f. for large 
values of D is given in Section 5.4 of Part II, describing the bending tests, which is 
to be published later. Fig. 5 also shows that the s.c.f. increases continuously as r/d 
is decreased. 


The curves shown as full lines in Fig. 6 have been obtained, by interpolation, 
from the combined results of Jacobsen’ and Wiegand’, which have been generally 
accepted as the design criteria for shouldered shafts subjected to torsion. These 
curves are in fair agreement with the results of the present tests (shown as chain- 
dotted lines) for high values of ¢) D, but show increasingly marked differences as 
d/D and r/d are decreased. 


Jacobsen made use of an electrical analogy to obtain results. for 0-005 < r/d 
< 0:12, which indicate that the s.c.f. is asymptotic to a limit as d/D decreases to 
zero, whereas the results of the present tests indicate that the factor reaches a 
maximum value and then decreases for d/D less than 0-7. The analogue requires 


*The author is indebted to Miss M. E. Michael, of the Technical Department of the Royal 
Aeronautical Society. for the preparation of these curves. 
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the use of a steel “razor blade” model (shown in Fig. 7) having the same profile as 
a longitudinal section through the shouldered shaft. It is necessary for the thick- 
ness of the steel model to vary as the cube of the radial distance from the shaft axis 
Since zero thickness is unobtainable in practice, Jacobsen used a model representing 
a 12 in. diameter shaft having a 14 in. diameter axial bore. The thickness of the 
model varied from 0-002 in. at the bore to | in. at the outer surface. Jacobsen 
mentions that it is difficult to maintain the profile of the cross section within the 
small limits required for the electrical tests. 


The analogue used gives the twisting moment as N = 2«/. where / is the total! 
current flowing through the model, and the shear stresses as 
ro ur CP / 


and zo=ure?/ 


where and ed /¢éz are the potential gradients in the appropriate directions 
on the surface of the model. The probes used for measuring the potential gradients 
were 0-1 in. apart (cf. 0-5 in. maximum shoulder root radius) and it was estimated 
that an error of up to 2 per cent was incurred by faulty probe contact. The 
potential difference across the probes was measured by a galvanometer, which had 
to be recalibrated immediately after each reading in order to obviate the effects of 
drift. The difference between the potential gradients measured at corresponding 
points on the curved and on the flat sides of the model was as much as 10 per cent 
when r/d or D/d became large. Jacobsen suggests that his results are subject to 


r/d 


0:05 


d/p 


Ficure 6. Comparison of the combined results of Jacobsen and Wiegand with those 
of the present investigation. 
Jacobsen®? (electrical analogue) and 
Wiegand? (precision strain gauge) 
-— — Allison. 
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FicurRe 7. Diagram showing the model used by Jacobsen, and the dimensions 
of the shaft which it represents. 


the greatest inaccuracy in these regions and it is interesting to note that the 
maximum divergence between his results and those of the present investigation 
occurs for small d/D ratios. 


Wiegand obtained results in the range 0°1 << r/ D < 0:25, using a strain gauge 
of 1-3 mm. gauge length. For large values of d/D,. the results of these tests are in 
agreement with those of the present investigation. Despite the fact that Wiegand’s 
experiments cover the range 0 << d/D <0°5, no decrease in the s.c.f. is predicted 
for the smaller values of d/D. From a sample calculation it appears that the 
maximum strain measured in the tests was of the order of 4:3« 10° (i.e. an 
extension of 2x 10 ° inches). It seems unlikely that the measurement of quantities 
of this order would lead to the calculation of stress concentration factors having an 
accuracy which is better than +10 per cent, although Wiegand estimates the 
accuracy of his curves to be +5 per cent. The difference in the results may be 
attributed to the fact that Wiegand used the base length of a mechanical extenso- 
meter to measure strains which actually occur on a curved arc. This would explain 
the divergence of the families of curves as r/d decreases and the fact that the 
position of the peak stress recorded by Wiegand is farther up the fillet than has been 
found in the present investigation. 


Undoubtedly the s.c.f. for small d/D ratios are too high to be reconciled with 
the photoelastic results, even if the maximum error band is put on the families of 
curves. and it can only be concluded that the results obtained in this region by 
Jacobsen and Wiegand are incorrect. 


It can be shown, in the particular case of the torsion of a cylindrical shaft of 
variable cross section, that the stress system is non-dilational and is therefore 
independent of the elastic constants. The implication of this result is that the stress 
systems in geometrically similar shafts, made of different materials, will be exactly 
comparable and that the stress concentration factors obtained using a photoelastic 
resin as the shaft material will be identical with those obtained in a similar 
steel shaft. 
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The only attempt to obtain a theoretical solution for the shouldered shaft in 
torsion is that due to Sonntag? who has taken a polynomial solution of the torsion 
function # in an attempt to solve the problem for a number of grooves and shoulders. 
He does not fully satisfy the boundary conditions in any of the cases and, for the 
shouldered shaft in particular, considers only the first four terms of the solution, 
without first establishing the convergence of the series. No reliance can be placed 
upon the results obtained by Sonntag. 
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Review 


Jet Propulsion Engines. Volume XII of High Speed Aerodynamics and Jet Propulsion. 
Edited by O. E. Lancaster. Oxford University Press, London, 1959. 799 pp 


Illustrated. £7 


Most readers must surely have met earlier volumes of this admirable series, whose 
only and inescapable fault is that they are so packed with information that their price 
individually, let alone in tote, is bound to deter many who would wish to possess them 
So it may be said at the outset that the volume under review is well up to standard. As 
a collective work, it is probably best to deal with it by taking the various contributions 
separately 

A. Historical Development of Jet Propulsion, by F. J. Malina, R. C. Truax and 
A. D. Baxter. is an admirably concise survey of the history of the various forms of 
reaction engine (including rockets, a legitimate inclusion under the heading of the book 
which is pursued later), and has a valuable bibliography 

B. Basic Principles of Jet Propulsion, by Maurice Roy, is the expected masterly 
disentanglement of the subtle and often misunderstood intricacies of thermal and 
propulsive efficiencies for various engine configurations. 

C. The Turbojet Engine. by C. A. Meyer, is a selection of topics under this heading, 
covering performance and control, thrust augmentation and testing: the diagrammatic 
presentation of the matching problem on page 128 repays close study. 


May 1961 199 


if 
| 


REVIEW 


D. The Turboprop Engine, by |. H. Driggs and O. E. Lancaster, not only presents 
the highly original method of performance analysis which the authors have made familiar 
in other publications, but also takes note—not before time—of the fact that the turbo- 
prop has a “prop” as well as a “turbo”, and spends useful time discussing propeller 
performance. 

E. The Ramjet Engine, by D. D. Wyatt and B. T. Lundin, is, again, chiefly 
concerned with performance aspects, including the very important topic of variation of 
specific heats, but it also gives a survey of the characteristics of test facilities. 

F. Intermittent Jets, by J. V. Foa, discusses not only the pulsejet, but also a variety 
of less-familiar proposals for cyclically-operating air-breathing engines. 

G. The Liquid Propellant Rocket Engine, by M. Summerfield, is less analytical 
than some of the preceding sections and goes more fully into design considerations 
Paradoxically, it is nevertheless devoid of worthwhile photographs of current practice, 
unlike the earlier parts. 

H. Solid Propellant Rockets, by C. E. Bartley and M. M. Mills. places more 
emphasis on combustion problems and grain design, as is to be expected, but again 
dwells briefly on mechanical design. 

I. The Ramrocket, by 1. Glassman and J. V. Charyk, is necessarily a mainly 
theoretical survey. 

J. Jet Rotors, by A. Gail, discusses the dynamics and aerodynamics of jet-driven 
helicopter and propeller blades. 

K. Atomic Energy in Jet Propulsion, by R. Zirkind surveys, in turn, nuclear 
physics, reactor principles, shielding and cooling. 

L. Future Prospects of Jet Propulsion, by F. Zwicky, applies morphological and 
topological principles to the discussion of possible future forms of engine, and will 
certainly cause most readers to wonder what has hit them. 


It will have been seen that the mere statement of the various fields covered by this 
vast book has taken a lot of space. An overall assessment is very hard to make fairly 
One feels that considerations of space and time have militated against any close integra- 
tion of the various topics covered, some of which should not have received such full 
treatment in a balanced review of the jet propulsion field at large: this is particularly 
true of Sections F, I, J, and K. Contrariwise, some mention of such things as ion rockets, 
turbo-rockets and the various forms of turbo-fan would have had more general appeal 
than the high-flown language of much of Section L. 


Were the last of these sections to be replaced by a more conventional peep into the 
future, and the rest coalesced into a shorter précis of fringe topics, one would have had 
a golden, rather than a curate’s, egg. As it is the authors have dealt nobly with an 
enormous field, and their joint labours form a classic milestone. But what a pity the 
editor issued invitations to contribute so freely: and surely the title ought to be more 
honest about the primarily performance and theory bias of the whole work 


Nevertheless, no serious student of the propulsive art ought to neglect to dip into 
this book variously and often, for he will find much to inform, stimulate, puzzle and 
occasionally infuriate him 


J. R. PALMER. 
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The Bristol Siddeley Olympus turbojet powers the |The Bristol Siddeley Gamma rocket motor powers 
Avro Vulean V-bomber. the Saunders-Roe Black Knight space probe. 


The Bristol Siddeley Thor ramjet powers the Bristol/ The Bristol Siddeley Proteus turboprop powers the 


Ferranti Bloodhound guided missile. Bristol! Britannia military transport. 


} 


The Bristol Siddeley Marine Proteus powers the Bristol Siddeley Mavbach diesel engines power the 


“Brave” class 96-1t boats built by Vosper Ltd. British Railways “Bristolian” ex] 


The Bristol Siddeley BS 53 turbofan powers the The Bristol Siddeley Olympus will power BAC’s new 
Hawker P 1127 VTOL strike fighter. tactical /support reconnaissance aircraft—the TSR 2. 


Bristol Siddeley Engines Limited— 
one of the largest producers 
of motive power units in the worid 
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NOTICE TO CONTRIBUTORS 


1. All communications should be addressed to the Editor, The Aeronautical 
Quarterly, 4 Hamilton Place, London, W.1. 


Papers for consideration should be typewritten with double spacing and be 
accompanied by a summary of not more than 250-300 words. Papers must be 
as short as possible, otherwise considerable delay in publication will result, and 
in any case they should not be longer than twenty foolscap pages. Full 
descriptions and complete results of work, as written for an internal report, are 
not suitable for publication in The Aeronautical Quarterly, and such reports 
need to be rewritten in concise form. Only necessary equations should be 
included in the text; routine mathematics should be omitted or given as an 
Appendix to the paper. Illustrations should be reduced to a minimum. The 
titles of papers should be kept short—with sub-titles if necessary. 


3. Only simple symbols and formulae should be typewritten. All others should 
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Suffixes and prefixes must be indicated clearly and complicated suffixes 
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4. A list of all symbols used, both in the illustrations and text, must be given 
whether they are standard or not. 

5. Clear diagrams or sketches must be submitted for redrawing by the Society’s 
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Photographs must be not less than half plate size and must be clear black 
and white glossy prints. 
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